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I.  This  report  is  a final  report  covering  research  studies  during  the  five 
year  period,  I April  1972  to  30  March  1977,  on  extrinsic  silicon  infrared 
detectors.  The  main  objective  was  to  study  the  material  parameters  which 
limit  and  control  the  performance  of  silicon  extrinsic  infrared  detectors. 
These  parameters  include  carrier  mobilities  and  the  thermal,  optical,  and 
generation-recombination-trapping  rates  at  impurity  (sensing  levels) 
conpensation  and  defect  centers.  The  experimental  approach  made  use  of 
ary' extended  transient  Junction  current  and  capacitance  measurement  tech- 
niques to  determine  photo-ionization  cross  sections,  and  thermal  capture  and 
emission  rates  as  a function  of  temperature,  electric  field,  and  concentration 
of  impurity  and  defect  centers  as  well  as  to  determine  the  concentration  pro- 
file of  these  centers.  Theoretical  interpretation  included  calculations  of 
photo-ionization  cross  sections  and  ion  scattering  mobilities  using  the 
pseudo-impurity  potential.  A new  objective  added  during  the  fifth  year  of 
the  contract  became  the  central  force  of  the  follow-on  contract  (FI9628-77-C- 
0138);  namely,  the  study  of  the  physical  mechanism  and  circuit  model  of  low 
frequency  l/f  noise  generated  at  the  Si02-Si  interface  states  which  are 
relevant  to  signal  processing  applications.  Some  of  the  important  results 
include: 


a)  The  absence  of  impurity  ion  position  correlation  in  mobility  scattering. 

b)  Further  development  of  transient  capacitance  technique  so  as  to  delect 
and  characterize  minute  amount  of  impurities  and  defects  introduced  during 
high  temperature  processing  as  well  as  to  fine  probe  the  concentration  pro- 
files of  these  centers. 


c)  The  study  and  characterization  of  numerous  impurity  and  defect  centers 
in  silicon  including  phosphorous,  bismuth,  gold,  magnesium , silver,  cobalt, 
chromium,  sulfur,  tantalum,  zinc,  aluminum-oxygen  complexes,  thermally 
generated  donors,  and  defect  centers  generated  by  ion  implantation. 

Support  from  this  contract  was  acknowledged  in  6 master's  theses,  7 doctoral 
theses,  and  51  papers  published  in  the  open  literature  during  its  duration. 


2.  Silicon  extrinsic  detectors  are  finding  rapidly  widening  use  in  surveillance, 
guidance,  and  optical  communication  systems.  This  study  has  contributed  basic 
information  needed  in  extending  the  range  of  applicability  of  IR  detectors  through 
the  use  of  new  sensing  impurities.  In  addition,  it  has  provided  basic  material 
parameters  on  shallow  and  deep  level  impurity  and  defect  centers  which  control 
and  limit  the  sensitivity,  uniformity  and  bandwidth  of  presently  used  detectors 
and  projected  large  scale  arrays  based  on  these  detectors. 


ANDREW  C.  YANG 
Project  Engineer!^ 
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SUMMARY 


Studies  of  the  material  parameters  which  limit  the  performance  of 
silicon  extrinsic  infrared  detectors  are  reported.  These  parameters  include 
the  carrier  mobilities  and  the  thermal,  optical  and  Auger-impact  generation- 
recombination-trapping  rates  at  impurity  and  defect  centers  and  impurity-defect 
complexes.  The  transient  capacitance  methods  have  been  used  to  obtain  these 
material  characterization  parameters  as  well  as  the  generation-annealing- 
diffusion  kinetics  of  the  centers.  The  transient  capacitance  methods  have 
been  extended  and  refined  during  this  contract  to  allow  for  fine  probing  of  the 
concentration  profile  of  these  centers. 
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I. 


INTRODUCTION 


Infrared  detector  is  one  of  the  most  important  elements  in  optical 
communication  receivers.  Many  types  of  solid  state  detectors  (particularly 
semiconductor  detectors)  have  been  developed  to  meet  specific  system  re- 
quirements .[  I , 2 ] Both  the  junction  (photovoltaic,  PV)  and  bulk  (photoconduc- 
tive,  PC)  type  of  devices  have  been  widely  used.  The  p-n  junction  type  gen- 
erally operates  as  an  intrinsic  detector  (optical  generation  of  electron-hole 
pairs  across  the  band  gap) . The  Schottky  barrie"  junction  type  operates  by 
collecting  the  optically  excited  electrons  from  the  metal  Fermi  surface  over 
the  metal-semiconductor  barrier.  The  photoconductive  type  can  be  either  an 
intrinsic  detector,  employing  interband  optical  generation  of  electron-hole 
pairs  in  relatively  pure  semiconductor  bulk,  or  an  extrinsic  detector,  using  a 
doped  semiconductor  in  the  deionization  temperature  range  where  a fraction 
of  the  impurities  are  neutral  so  that  the  trapped  electrons  (or  hi.xes)  can 
be  optically  excited  in  the  detection  process.  Germanium,  doped  “.fith  Au, 

Zn,  Hg,  etc. , was  widely  used  as  cooled  extrinsic  photoconductive  detectors 
in  the  intermediate  infrared  range  in  the  past.  More  recently,  intrinsic 
photovoltaic  and  photoconductive  detectors  in  HgCdTe  have  been  developed. 

From  the  point  of  view  of  ultimate  detector  sensitivity  or  detectivity 
D*,  and  the  signal  to  noise  ratio,  there  appears  to  be  little  differences 
among  the  various  materials  when  the  detector  is  operated  under  background 
limited  conditions.  However,  in  order  to  reach  the  photon  background  noise 
limited  condition,  the  intrinsic  noise  sources  within  the  detector  must  be 
minimized.  It  is  felt  that  silicon  extrinsic  detectors  may  offer  some  advan- 
tages in  reducing  and  controlling  the  noise  sources  in  individual  detector 
devices  and  in  array  and  integrated  structure  applications  due  to  the  well 
developed  silicon  process  technology.  In  addition,  there  are  also  some  ad- 
vantages In  radiation  hardening  from  the  knowledge  accumulated  on  the  effect 
of  ionizing  radiation  and  radiation  defects  on  silicon  device  characteristics. 

This  study  is  concerned  with  the  theoretical  and  experimental  inves- 
tif'.ation  of  rivombinat  i on-gent^rat  i on-t  rapping  kinetics  and  the  fundamental 
properties  of  tlie  material  parameters  siich  as  the  thermal  capture  and  emission 
rates  of  trapped  electrons  and  holes  and  photoionizat ion  cross-sections  of 
impurity  centers  in  silicon  materials  useful  for  extrinsic  infrared  photodetector 
app] ications. 
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1.1  OBJIiCTIVnS 


The  performance  of  an  extrinsic  photoconduct ive  detector  Is  determined 
by  the  photoconductivity  and  the  detectivity  or  the  signal  to  noise  ratio. 

As  a first  example,  let  us  consider  an  n-type  silicon  without  deep  level  re- 
combination centers  whose  energy  level  diagram  is  given  in  Figure  1.1.1.  The 
photocurrent  is  given  by 


J = qvi  EAN 
S ^ n 


(1.1) 


where  p is  the  electron  mobility,  E is  the  bias  electric  field  and  AN  is  the 
n 

steady-state  photoelectron  concentration.  For  low  steady-state  level,  i.e., 

AN<<N,  , , 
dark 


AN  = e^N-t  = o^*!'  N„t  = a°4>  t = t /L 
nun  nsDn  nsn  sn 


(1.2) 


Here,  is  the  photon  flux  (photons/cm^-sec) , o°  is  the  average  photoioniza- 
tion cross  section  over  the  photon  spectrum  and  the  detector  spectral  range, 
a°  = is  the  average  absorption  coefficient  (1/cm),  n = a°L  is  the  average 

internal  quantum  efficiency,  L is  the  sample  length  In  the  direction  of  the 

illumination,  and  N_  = N„..  “ N.  . is  the  dark  trapped  electron  concentration  at 
D DD  AA 

the  donors. 

The  low  level  steady-state  electron  lifetime  is  given  by 


T 

n 


(1.3) 


where 


P = N - N 
D DD  D 


AA 


(l.A) 

(l.Aa) 


For  the  ideal  example,  the  slgnal-to-noise  ratio  and  the  detectivity 
are  given  by 


/dIS/N  = ())«>  /2)/n„'  T A/P,^N_L 
s DD  n D D 


(1.5) 


and 


D*  . (,„„/2.,.^)4c‘/I.)(1  - 1/c, 


-2- 


(1.6) 


Figure  1.1.1  Eaergy  band  diagram  of  a n-type  semiconductor  doped  with 
a shallow  level  acceptor  with  energy  level  Ep^  and  shallow 
level  donor  witn  energy  level  Ed.  The  thermal  and  optical 
transition  rates  are  indicated  next  to  the  arrows. 


Generated  by  ionizing  rodiation^ 


Figure  1.1.2  Energy  level  diagram  of  a n-type  semiconductor  doped  with 
a shallow  level  acceptor  with  energy  level  E;^,  a shallow 
level  donor  with  energy  level  Eq  and  a deep  recombination 
impurity,  such  as  gold,  with  an  energy  level  of  E/^y.  The 
transition  rates  are  indicated  next  to  arrows. 


Here  Af  is  the  bandwidth,  A is  the  active  detector  area,  hv^  is  the  photon 
energy  and  c^  is  the  thermal  capture  rate  (cm^/sec)  of  electrons.  In  addi- 
tion, the  3 d.b.  bandwidth  is  given  by 


tDT  = 1 
n 


where  u is  the  signal  frequency. 

Based  on  these  results  for  an  ideal  detector  with  no  deep-level  re- 
combination centers,  it  is  evident  that  the  material  parameters  which  limit 
the  performance  of  n-type  extrinsic  detectors  are: 

(i)  Mobility  of  marjority  carriers,  [see  Equation  (1.1)]. 

(ii)  Photoionization  cross-section  of  the  majority  carriers,  o°  [see 
Equation  (1.2)]. 

(iii)  Thermal  capture  rate  of  the  majority  carriers  at  the  shallow 
level  donor,  c^  [see  Equations  (1.3),  (1.5),  and  (1.6)]. 

In  practice,  deep  level  recombination  centers  are  always  present  in 
the  detector,  either  from  those  present  in  the  original  crystal  or  Introduced 
during  high  temperature  device  fabrication  processes.  The  presence  of  these 
deep  level  centers  drastically  affects  the  detector  sensitivity  and  bandwidth. 
As  an  example,  let  us  consider  the  presence  of  the  gold  acceptor  recombina- 
tion center  in  n-Si  whose  energy  level  diagram  is  given  in  Figure  1.1.2.  The 
parameters  that  determine  the  photoelectron  concentration  in  Equation  (1.2) 
are  now  given  by 


«D  = 


DD 


Au 


(1.7) 


and 


T 

n 


l/c^^P, 
n D 


(1.8) 


= 1/c  (N. . + N.  ) 
n AA  Au 


(1.9) 


Thus,  the  photoconductivity,  which  is  proportional  to  N^j  as  indicated  in 
Equation  (1.2),  is  reduced  by  the  compensating  gold  acceptor  centers  as  indi- 
cated in  Equation  (1.7).  In  addition,  the  bandwidth  is  increased  by  both  the 
shorter  electron  lifetime  as  indicated  in  Equation  (1.9)  and  the  additional 

recombination  of  electron  and  hole  through  the  gold  recombination  centers  at 

t -1  t t 

a rate  of  c,  N.  , giving  a total  reciprocal  lifetime  of  t = C (N.,+N,  )+C.  N.  . 
^ Au  Au’  06  r ^ ^ n AA  Au  Au  Au 

Here  C.  is  signal  level  dependent  and  reaches  C.  C.  / (C . ) at  high 
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Based  on  the  analyses  of  the  two  examples  just  given,  the  objectives 
of  this  contract  were  delineated.  They  are  summarized  below: 

(1)  To  study  the  effect  of  majority  and  minority  carrier  lifetimes 
and  mobilities  in  the  presence  of  recombination  centers  such  as  gold,  radia- 
tion damage,  thermally  generated  defects  and  others; 

(2)  To  make  uje  of  and  to  extend  the  transient  junction  current  and 
capacitance  measurement  techniques  to  determine  the  photoionization  cross 
section,  and  thermal  capture  and  emission  rates  as  a function  of  temperature, 
electric  field  and  impurity  concentration  at  shallow  and  deep  level  impurity 
centers,  as  well  as  to  determine  the  concentration  profiles  of  these  centers; 
and 

(3)  To  perform  theoretical  calculations  and  interpretations  of  the 
photoionization  cross  sections  and  ion  scattering  mobilities  using  the  pseudo- 
impurity  potential  models. 

In  addition  to  these  original  objectives,  a new  objective  was  added 
during  the  fifth  year  on  the  study  of  the  physical  mechanisms  and  circuit 
models  of  low-frequency  1/f  noise  generated  at  the  Si02-Si  interface  states 
which  are  relevant  to  MOS  transistors  and  CCD  devices  in  image  processing 
applications. 

1.2  CONTENT 

This  final  technical  report  will  give  a summary  of  those  projects 
which  have  been  reported  as  journal  articles  in  open  literature.  A detailed 
account  of  unpublished  results  will  also  be  given. 

The  work  on  electron  and  hole  mobilities  in  silicon  is  given  in 
chapter  II.  The  following  format  of  subsections  will  be  followed  for  this 
and  subsequent  chapters.  The  problem  statement  is  given  in  section  2.0.  The 
experimental  techniques  are  described  in  section  2.1.  The  experimental  results 
are  given  in  the  subsections  of  section  2.2.  The  electron  mobility  is  given 
in  section  2.2.1  and  the  hole  mobility  in  section  2.2.2.  Theoretical  analysis 
and  interpretation  are  given  in  section  2.3.  Summary  and  suggestions  for  fu- 
ture work  are  given  in  section  2.4. 

The  work  on  the  measurements  and  interpretation  of  the  properties  of 
recombination  centers  is  covered  in  chapter  III.  Concluding  remarks  are  made 
in  chapter  IV. 


1.3  STATISTICS 

During  the  five-year  contract  period  (01APR72  to  31MAR77),  twelve 
(12)  post-doctoral  associates  were  partially  supported  by  the  contract  and  have 
contributed  to  the  technical  results  reported  here.  The  names  of  these  asso- 
ciates and  their  current  affiliations  are  given  in  Appendix  A. 

Twenty-one  (21)  graduate  students  were  partially  supported  by  this 
contract,  resulting  in  six  (6)  master  theses  in  electrical  engineering,  four 
(4)  doctoral  theses  in  electrical  engineering,  four  (4)  doctoral  theses  in  physics, 
and  two  additional  doctoral  theses  near  completion.  The  names  of  these  students, 
their  thesis  titles,  degree  dates  and  current  affiliations  are  given  in  Appendix 
B. 

There  were  a total  of  fifty-one  (51)  papers  published  in  the  open 
literature  during  the  duration  of  this  contract  and  one  (1)  before  the  date  of 
this  report.  Tnese  are  listed  in  Appendix  C. 

II.  ELECTRON  AND  HOLE  MOBILITIES 


2.0  PROBLEM  STATEMENT 

We  have  already  seen  in  section  1.1  that  photoconductivity  is  directly 
proportional  to  the  mobility  of  the  majority  carriers  at  low  level.  It  is  pro- 
portional to  the  sum  of  the  mobilities  of  the  majority  and  minority  carriers  at 
high  level.  Thus,  higher  photoconductivity  and  signal  would  result  when  mobil- 
ities are  higher.  However,  the  detector  sensitivity  is  higher  if  the  concentra- 
tion of  the  optically  active  impurity  center  is  higher.  But  large  impurity 
concentration  will  reduce  carrier  mobility  due  to  Coulomb  scattering  by  the 
randomly  located  impurity  ions.  Thus,  there  is  an  optimum  impurity  concentration 
for  maximum  detector  sensitivity.  This  optimum  concentration  may  be  further 
affected  by  three  factors:  (1)  Maximum  solubility  of  the  Impurity,  (2)  Posi- 
tion correlation  of  impurity  atoms,  and  (3)  Mobility  in  impurity  bands.  The 
first  factor  may  make  the  impurity  concentration  at  maximum  sensitivity  unreach- 
able. The  second  factor  will  increase  the  sensitivity  but  is  limited  by  funda- 
mental thermodynamic  and  statistical  fluctuations  during  the  high  temperature 
fabrication  processes.  The  third  factor  is  important  when  the  impurity  concen- 
tration is  higher  than  about  10^®  atoms/cm®.  This  range  is  important  in  certain 
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applications  (e.g.,  emitter  or  bipolar  transistor)  but  is  unlikely  to  be  very 
important  in  extrinsic  detector  unless  substantial  position  correlation  can  be 
maintained  to  avoid  impurity  band  when  the  concentrations  of  optically  sensitive 
impurity  centers  are  higher  than  about  10^^  atoms/cm^. 

The  purpose  of  this  study  is  to  develop  a theory  that  can  take  into 
account  the  position  correlation  of  the  scattering  centers  and  to  determine  if 
there  is  any  position  correlation  in  samples  prepared  by  conventional  techniques. 
A new  mobility  measurement  technique  is  also  developed  in  this  study. 

2.1  EXPERIMENTAL  TECHNIQUES 

The  conventional  mobility  determination  technique  gives  Hall  mobility. 
It  is  not  possible  to  convert  Hall  mobilities  into  conductivity  mobilities  with 
a certainty  better  than  a factor  of  25%  and  often  as  much  as  2 due  to  the  un- 
known temperature  and  energy  dependences  of  the  scattering  mechanisms.  The 
problem  lies  in  the  measurement  of  the  carrier  concentration,  N,  when  calculating 
the  mobility  from  the  measured  conductivity,  a. 


Hall  coefficient  measurement  gives  a carrier  concentration  which  is  subject  to 
the  large  uncertainty  discussed  above.  An  alternative  technique  is  developed 
to  measure  the  carrier  concentration,  using  the  capacitance-voltage  character- 
istics of  a reverse  biased  Schottky  barrier  diode. [3]  However,  this  method  is 
subject  to  some  uncertainty  at  high  impurity  ion  concentrations  and  at  low  tem- 
peratures due  to  impurity  deionization  effect. [4]  What  is  measured  from  the 
C-V  curves  is  the  total  net  impurity  concentration  (donor  minus  acceptor)  and 
not  the  carrier  (electron)  concentration.  Unless  all  the  impurities  are  ionized, 
N = Npjj  - does  not  hold  but  N < so  that  the  calculated  mobility 

would  be  higher  than  the  true  mobility. 

2.2  EXPERIMENTAL  RESULTS 

The  C-V  technique  is  employed  to  obtain  the  hole  mobility  in  p-type 
silicon  and  this  is  reported  in  section  2.2.1.  Experimental  data  of  electron 
mobility  are  obtained  from  the  open  literature  which  are  mostly  Hall  mobility 
or  effective  conductivity  mobility  without  the  correction  for  the  deionization 
and  impurity  band  effects. 


t 


2.2.1  HOLE  CONDUCTIVITY  IN  p-Si 

The  results  were  reported  In  two  published  papers.  The  experimental 
technique  and  the  detailed  results  were  given  by  Tsao  and  Sah. [3]  The  analyses 
made  by  Woodley  and  Sah,  Including  the  deionization  effect,  are  summarized  in 
Table  1 and  Figures  2.2.1,  2.2.2,  2.2.3,  and  2.2.4. [4]  The  effect  of  deioniza- 
tion on  the  hole  concentration  calculated  from  the  C-V  data  Is  shown  In  Figure 
2.2.1  where  the  dots  are  the  net  Impurity  concentration  and  the  solid  curves 
are  the  hole  concentrations  for  i.  -’s  p-Sl  samples  Investigated  whose  net  Ionized 
Impurity  concentration  Is  shown  ’able  1.  The  corrected  (solid  curves) 
and  uncorrected  (dots)  mobilities  'or  two  samples  are  shown  In  Figure  2.2.2  as 
a function  of  temperature.  The  corrected  mobilities  with  deionization  taken 
Into  account  are  shown  in  Figure  2.2.3  for  all  six  samples  where  the  lattice 
scattering  mobility  Is  also  shown  by  the  dashed  curve.  A comparison  of  the 
room-temperature  hole  conductivity  mobility  as  a function  of  boron  Impurity  con- 
centration obtained  In  this  work  with  those  of  Irvin  [5]  Is  given  In  Figure 
2.2.4.  The  present  data,  after  making  the  deionization  correction,  compares 
favorably  with  our  previously  unpublished  data  while  It  deviates  significantly 
from  the  frequently  used  Irvin  curve. 


2.2.2 


ELECTRON  MOBILITY  IN  n-Sl 


The  data  of  the  electron  mobility  were  collected  from  open  literature. 
Most  of  these  are  Hall  mobilities,  e.g. , the  electron  concentration  Is  calculated 
from  the  Hall  coefficient  data.  Some  of  the  mobilities  at  very  high  Impurity 
concentrations  are  obtained  by  assuming  that  the  electron  concentration  is  equal 
to  the  majority  impurity  (phosphorus  or  arsenic)  concentration.  The  mobility 
due  to  Ion  scattering  alone  was  then  computed  from  the  Matthlessen  rule 

^^^I  ^^’’'Experimental  ^ ^^^L 

where  the  lattice  scattering  mobility,  )Jj^  = 1330  cm^/V-S  at  300‘’K,  is  the  exper- 
imental mobility  asymptotically  extrapolated  to  low  ion  concentration.  These 
experimental  results  are  given  in  Figure  2.2.5  at  300K  and  Figure  2.2.6  at  77K. 

2.3  THEORETICAL  ANALYSIS  AND  INTERPRETATION 

The  conventional  theories  of  ion  scattering  (Brook-Herring  and  Conwell- 
Welsskoff  formulae)  are  based  on  the  model  of  a completely  random  distribution 


-8- 


Table  1 


Figure  2.2 


Impurity  Concentration  of  p-type  Silicon 


Device 

^AA- 

68C24 

9.71 

X 

10^^ 

7412 

5.52 

X 

10^5 

7645 

1.06 

X 

10^^ 

7735 

2.27 

X 

10^^ 

7914 

2.79 

X 

10^" 

8053 

4.71 

X 

10^^ 

TE»/P£RATURE  (K) 


1 Ionized  impurity  and  hole  concentrations  versus  temperature 
for  devices  used  in  [3].  Dots  are  values  of  impurity  concen- 
tration obtained  from  C-V  measurements.  Solid  curves  are  hole 
concentration  corrected  for  deionization  effects. 


g id* 
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TEMPERATURE  (K) 


Figure  2.2.2 


Conductivity  mobility  of  holes  versus  temperature.  Points 
represent  uncorrected  values  for  lowest  (10l5B/cm3,  dots)  and 
highest  (IxlO^^B/cm^,  circles)  impurity  concentrations  used  in 
[3].  Solid  curves  are  corrected  for  deionization  effects.  Dashed 
curve  represents  mobility  in  the  absence  of  impurity  scattering. 
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Figure  2.2.3  Conductivity  mobility  of  holes  versus  temperature  for  devices 
used  in  [3]  corrected  for  deionization  effects.  Dashed  curve 
represents  mobility  in  the  absence  of  impurity  scattering. 


Figure  2,2.4  Conductivity  mobility  of  holes  versus  Impurity  concentration  at 

room  temperature.  Results  of  this  work  (points)  compared  to  values 
calculated  from  fits  to  published  resistivity  versus  Impurity 
concentration  data  by  Irvin  [5]  (solid  curve)  and  previous  data 
used  by  Sah  (dashed  curve).  The  estimated  accuracy  Is  about  15%. 


fji.  Mobility  (cmVv-sec) 


Figure  2.2.5  Theoretical  ionized  impurity  scattering  electron  mobilities  in  silicon 
at  300°K  as  a function  of  ion  concentration.  Pgjj  is  the  Brook-Herring 
theory  without  any  ion  position  correlation  and  IIq  is  the  Brook-Herring 
theory  modified  to  include  ion  position  correlation.  Debye-Huckel 
screening  at  800°K  is  assumed  for  ion  position  correlation. 


Figure  2.2.6  Theoretical  ionized  impurity  scattering  electron  mobility  in  silicon 
at  77°K  as  a function  of  ion  concentration.  Hard  sphe.re  ion  position 
correlation  is  assumed  and  electron  screening  given  by  equations 
(2.9),  (2.11)  and  (2.12)  are  used. 


of  point  ions.  This  model  neglects  any  position  correlation  of  the  ions.  Po- 
sition correlation  gives  rise  to  coherency  and  hence  less  random  scattering  or 
higher  mobility. 

An  isotropic  elastic  scattering  theory  is  worked  out  which  takes  into 
account  position  correlation  of  the  Coulomb  scatterers.  The  following  assump- 
tions are  made:  Born  approximation,  linear  dielectric  and  Debye-Huckel  screen- 
ings, and  hard  sphere  position  correlation  with  radius  r^  given  by  the  average 
concentration  of  the  scatterer  or  the  Debye  screening  length  for  position 
correlation. 

The  theory  does  not  make  an  assumption  about  the  electron  concentra- 
tion so  that  it  applies  to  both  the  nondegenerate,  intermediate,  and  degenerate 
ranges.  However,  in  the  high  concentration  degenerate  or  near  degenerate  range, 
the  concentration  of  the  scatterers  is  so  high  that  some  of  them  are  not  ionized. 
In  addition,  the  semiconductor  conduction  band  edge  is  no  longer  well  defined 
and  impurity  band  conduction  may  not  be  negligible.  These  factors  are  not 
taken  into  account  so  that  the  theory  is  not  reliable  in  the  high  concentration 
range. 

The  result  of  the  theory  is  given  below  and  a family  of  curves  is  nu- 
merically computed  using  an  IBM  360  machine  and  presented  in  Figures  2.2.5  and 
2.2.6.  These  are  briefly  discussed  below  and  compared  with  the  experimental 
data  described  in  section  2.2.2.  MKS  units  are  used  except  that  length  is  in 
cm  in  the  final  numerical  results. 


The  mobility  is  given  by 


Texp(X-X„)X^^^dX 
r 

[l+exp(X-Xp)]2 


f 

(l+exp(X-X.)J 
0 ^ 


(2.1) 


and  the  carrier  concentration  is  given  by 


N - »c’'l/2<V 


(2.2) 


Here,  e is  the  magnitude  of  the  electron  charge,  m^  is  the  conductivity 

effective  mass,  X • E/kgT  is  the  normalized  electron  kinetic  energy  with  the 

conduction  band  edge  = 0 as  the  reference,  E = Hxti,  where  m,  is  the 

L.  da 

density  of  state  effective  mass,  X = E^/k  T is  the  normalized  Fermi  energy, 

3/2  ^ r D 

Nc  = (g^/Aii^)  (2inn^kgT/h2)  is  the  effec  tive  density  of  state,  g^  is  the 
number  of  conduction  band  valley  and  Eg(Xp)  is  the  Fermi-Dirac  integral  of  the 
S-th  order  defined  by 


i 
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[l+exp(X-Xp)] 


(2.3) 


whose  derivative  is 


F' <X^)  - dFj/dXp  - F3_^(Xj) 


(2.4) 


The  reciprocal  momentum  relaxation  time  due  to  isotropic  elastic 
Coulomb  scattering  of  electrons  by  a partially  random  (finite  position  corre- 
lation) ion  charge  density  inhomogeneity  n^(X)  is  given  by 


8ire‘*m  k 
d 


(I~cos0)sin0d0 

(q)q‘* 


< nj^(q)n*(q)U  > 


(2.5) 


Here,  k is  the  wave  vector  of  the  incident  electron  wave  and  the  scatterer 
electron  wave  vector  is  k+q.  Since  the  scattering  is  elastic,  |k|  = jk+ql  or 
q =»  2k  sin(0/2)  where  0 is  the  scattering  angle  /k+q ,k  . Planewave  or  orthogo- 
nalized  planewave  representation  of  the  electron  wavefunction  is  used. 

The  dielectric  function  e(q)  contains  two  parts:  e(q)  = e^(q)  + ^^(*1) 
where  e-^(q)  comes  from  dielectric  screening  of  the  ion  charges  by  the  valence 
band  electrons  and  e^(q)  from  Debye  screening  by  the  conduction  band  electrons. 

e^(q)  is  nearly  given  by  the  static  dielectric  constant  [e^(0)  = = 

11.7  X 8.85  X 10~^**  F/cm  for  silicon]  since  the  dominant  scattering  events  off 
a Coulomb  scatterer  are  small  angle  scattering  0=0  and  in  nondegenerate  or  near 
degenerate  semiconductors,  k<<2v/a  where  'a'  is  the  lattice  constant.  If  it  is 
necessary  to  consider  laige  q events*  e^(q)  can  be  readily  obtained  from  band 
theory  calculations  which  are  now  available.  It  has  the  approximate  form  of 


e^(q)  = l+(E^-l)[l+(q2+2q-G-2q-k+2k-G)/G2]"2 


(2.6) 


where  G is  a reciprocal  vector  of  the  lattice.  This  shows  that  e (q)-<^:  at 

V O 

q-K)  or  q<<G  and  e^(q)-^l  when  q>>G. 

The  Debye  screening  function,  e^(q),  from  conduction  band  electron 
has  an  explicit  form  in  the  linear  screening  theory  given  by  [6] 


(2.7) 


where  the  square  of  the  Debye-Huckel-Thomas-Fermi  screening  wave  number,  q^, 
or  the  reciprocal  Debye  screening  length  = q^  is  given  by  (7,8) 
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(2.9) 


(e^m./ire  ) (3N/it) 
a o 


(2.10) 


If  large  q scattering  cannot  be  neglected,  an  explicit  form  of  is 


given  by  [6] 


o o I 2k^+q 

Ijj  = [e  N/(2Ep/3]{-j  + gj^q  *’"l  2ky-ql  ^ 


(2.11) 


where  and  = (3Tr^N/g^)^'  . 

One  of  the  problems  associated  with  the  theory  of  screening  of  the 
ion  charge  by  conduction  electrons  is  the  number  of  electrons  in  a Debye  sphere 


of  radius  q~'  or  Debye  screening  length,  » 9^^*  Evidently,  the  screening 


concept  breaks  down  if  there  are  few  electrons  in  a Debye  sphere.  Thus,  an 


alternative  form  for  q^  or  X^  is  to  take  it  as  the  average  inter-conduction -band- 


electron  spacing 


h “ •’d^  “ (3/47tN)^ 


(2.12) 


This  would  be  the  lower  limit  of  X^  below  which  the  simple  screening  model  is 


no  longer  meaningful. 


In  (2.5)  U is  the  volume  of  the  crystal.  n^(q)  is  the  Fourier  trans- 


form of  n^(r)  which  is  given  by  n^(r)-Nj  where  nj(r)  is  the  local  density  of 
the  ion  scatterers  while  =<nj(r)>is  the  average  ion  density.  < > is  the 

ensemble  average  over  all  configurations  of  ion  density  distribution.  < > 

could  also  be  evaluated  in  the  original  form  of  a double  sum  of  expflq* (R^-R^, ) ) 


over  all  R^  ion  coordinates.  For  the  hard  sphere  model  in  which  it  is  assumed 


that  no  two  ions  can  be  closer  than  a distance  of  r , then 

o 


<n^(q)n*(q)U>  = Nj,[l-3U^Nj(sinqr^-qr^cosqr^) / (qr^)  (2.13) 


where  is  the  hard  sphere  volume  given  by  = 4iir^/3.  The  product  “ C 

is  the  /ton  poi^iX^on  COfiAttcution  ^acXon.  in  the  hard  sphere  model.  When  C=l,  we 
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r = (3/4nN  )■ 
o 1 
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is  the  closest  approach  in  the  hard  sphere  model.  The  appropriate  value  of  C 
for  a particular  sample  depends  on  the  condition  (such  as  temperature  and  time 
of  diffusion  and  drift)  under  which  the  impurity  ions  are  introduced  during 
the  high  temperature  fabrication  processes.  C is  only  a qualitative  measure 
since  in  an  actual  case,  the.  exact  distribution  of  the  inhomogeneity  probably 
differs  from  the  hard  sphere  model. 

To  get  a more  realistic  model,  the  linear  screening  theory  is  also 
applied  to  the  ion  distribution  at  the  oxidation,  diffusion  or  crystal  growth 
temperature,  T^,  at  which  the  ions  are  introduced.  For  this  case, 

<n.(q)n*(q)„>  = Nj{l  - q^j/lq^  + 2:^0 ’ ^ 

where  q^,  * (e^N^/e  h T ) and  q^  = e^N  /c  k_T  . N is  the  average  concentra- 
ol  I o B o oa  a o B o a ® 

tion  of  the  a-th  species  (a  = donor  ion,  acceptor  ion,  electrons  and  holes)  and 
is  the  static  dielectric  constant  at  the  impurity  ion  introduction  tempera- 
ture T . 
o 

Numerical  calculations  are  carried  out  using  an  IBM  360  machine.  C 


is  taken  as  the  constant  parameter  and  is  assumed  equal  to  N in  (2.11)  in 

order  to  get  a r^  which  is  valid  only  if  all  the  impurity  atoms  are  ionized. 

E;^(q)  = E^(0)  = Eq  Is  assumed  in  (2.6)  which  is  a good  approximation  since  only 

small  q scattering  events  are  important. 

Families  of  curves  are  computed  for  two  temperatures  to  illustrate  the 

magnitude  of  the  correlation  effect.  For  mobility  at  300 ®K.  given  in  Figure  2.2.5, 

the  nondegenerate  Debye-Huckel  screening  given  by  (2.9)  is  used  for  both  e^(q)  in  (2.7) 

and  ion  position  correlation  given  by  (2.14).  A "fabrication"  temperature  of  T = 800“K 

is  assumed  where  p =n  =10  cm.  m =m  and  m.  = 0.323m  are  used.  The 
*^0  0 c d 

results  are  shown  in  Figure  2.2.5.  p is  the  Brook-Herring  theory  without  ion 

dH 

position  correlation,  i.e.,  <n.(q)n*u^N  , and  u_  is  that  with  correlation  given 

X X.  X (w 

by  (2.14).  The  nondegenerate  approximation  breaks  down  at  an  electron  concen- 
tration of  N = N„  = Np  = 2.8  X 10^^  cm“^.  The  ion  scattering  mobility  in  Figure  2.2.5 

—1  —1  2 

are  calculated  from  p_=  (p,  - p.  ) where  p^=  1330  cm  /V-S  and  p,  are  taken  from  (i) 

I data  L L data 

recent  Hall  mobility  data  of  Omelyanovskii,  Fistful  and  Milvidskli  [9],  (ii)  Irvin's 
resistivity  curve  [5],  and  (iii)  a conductivity  v.s.  impurity  concentration 
curve  obtained  by  Sah  in  1959  which  was  fitted  to  all  the  published  data  up  to 
that  date  [ 10] . 

Two  points  may  be  noted:  (1)  The  two  theoretical  curves  show  little 
effect  of  correlation  at  impurity  concentration  less  than  about  lO^^cm”^  and 


correlation  gives  only  about  1.5  times  higher  mobility  (u„  = 1.5m_„)  at 

C Bn 

Npp  = 2 X lO^^cm”^,  and  (2)  the  experimental  data  slightly  favors  the  Brook- 
Herring  model,  Mgy,  which  contains  zero  correlation.  The  latter  is  expected 
since  the  impurity  ions  are  introduced  into  these  samples  during  crystal 
growth  (-1420°C)  where  the  electron  concentration  is  so  high  that  all  long 
range  correlation  effects  are  screened  out. 

Theoretical  calculations  are  also  made  at  77K  and  are  shown  in 
Figure  2.2.6  and  compare  with  some  experimental  data  in  open  literature. [9,11] 
The  hard  sphere  model,  given  by  Equation  (2.13),  for  ion  position  correlation 
is  used.  Three  screenimg  models  are  used  in  the  theory  which  are  given  by 
Equations  (2.9),  (2.11),  and  (2.12). 

The  theoretical  curves  of  C » 0,  0.2,  0.4,  0.6,  0.8,  and  1.0  labeled 
(2.9)  in  Figure  2.2.6  use  the  nondegenerate  Debye-Huckel  screening  model  given 
by  Equation  (2.9). 

The  single  theoretical  curve  labeled  (2.11)  in  Figure  2.2.6  uses  the 
degenerate  Thomas-Fermi  screening  model  given  by  Equation  (2.11).  This  is  in- 
cluded to  illustrate  the  effect  of  large  q scattering. 

The  third  family  of  curves  with  C = 0,  0.2,  0.4,  0.6,  0.8,  and  1.0 
and  labeled  (2.12)  in  Figure  2.2.6  uses  the  average  interconduction  electron 
spacing,  given  by  Equation  (2.12)  as  the  screening  length.  It  is  Included 
again  only  as  an  illustration  since  this  model  does  not  Include  the  long  range 
property  of  the  Coulomb  force  in  the  electron-electron  interaction  and  hence 
underestimates  the  screening  effect,  resulting  in  lower  mobilities. 

The  comparison  of  the  theoretical  results  given  by  Equation  (2.9)  with 
the  experimental  data  of  Omelyanovskii  [9]  and  Chapman,  Tuft,  Zook  and  Long  [11] 
again  favor  the  zero  position  correlation  model  (C»0). 

The  theory,  given  by  curves  labeled  (2.9)  in  Figure  2.2.6,  does  show 
that  if  substantial  ion  position  correlation  can  be  achieved,  such  as  C-1,  then 
the  mobility  can  be  increased  by  a factor  of  four  over  that  of  a completely 
random  ion  distribution  at  an  ion  scatterer  concentration  of  =■  10^^  ions/cm^ . 

2.4  SUMMARY  AND  SUGGESTIONS 

A new  conductivity  mobility  measurement  technique  was  introduced  in 
which  the  carrier  concentration  is  determined  from  capacitance-voltage  measure- 
ments on  Schottky  barrier  diodes.  This  method  does  not  have  the  uncertainty 
inherent  in  the  Hall  effect  measurement  of  the  carrier  concentration  due  to  the 
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unknown  scattering  properties.  However,  the  C-V  curve  gives  the  net  ionized 
impurity  concentration  rather  than  the  carrier  concentration.  Thus,  a theo- 
retical correction  of  the  amount  of  impurity  ion  deionization  must  be  made. 

This  correction  is  negligible  at  room  temperature  and  low  impurity  concentra- 
tion but  becomes  appreciable  at  low  temperatures  and  high  impurity  concentra- 
tions. 

The  hole  mobility  data  obtained  from  the  C-V  technique  shows  that  the 

currently  accepted  values  due  to  Irvin  are  significantly  in  error  at  impurity 

d 5 3 

concentrations  above  10  Boron/cm  at  room  temperature  where  a small  impurity 
deionization  effect  can  be  corrected  accurately  in  our  data. 

Additional  experimental  work  using  the  Schottky  barrier  C-V  tech- 
niques appears  desirable  for  p-type  as  well  as  for  n-type  silicon,  the  latter 
not  studied  during  this  contract,  in  order  to  obtain  a reliable  set  of  con- 
ductivity mobility  versus  impurity  concentration  curves  over  the  practical  tem- 
perature range  of  -55°C  to  125°C  and  the  theoretically  interesting  range  of 
4®K  to  AQO^K.  These  data  would  provide  further  insight  into  the  scattering 
model  and  the  possibility  of  position  correlation  effects  which  would  give  high 
mobility  due  to  reduced  random  ion  scattering  as  well  as  high  detector  sensi- 
tivity. The  more  reliable  data  would  also  stimulate  additional  detailed  theo- 
retical work  in  the  mathematical  description  of  the  scattering  models  using 
impurity  pseudo-potential  theory. 

Further  studies  also  appear  desirable  to  determine  the  experimental 
methods  and  conditions  to  give  large  ion  position  correlation.  Large  ion  posi- 
tion correlation  would  reduce  scattering  and  hence  permit  the  use  of  materials 
with  high  concentrations  of  optically  sensitive  impurity  ions  to  achieve  high 
infrared  detector  sensitivity  as  well  as  high  mobility  which  further  improves 
photo  signal. 


III.  PROPERTIES  OF  IflPERFECTION  CENTERS  IN  SILICON 

3.0  PROBI.KM  STATEMENT 

It  was  shown  in  section  1.1  that  the  fundamental  parameters  which 
characterize  optically  active  impurity  or  defect  centers  for  detector  appli- 
cations are  the  thermal  capture  and  emission  rates  of  electrons  and  hole,  c^, 

t t t ^ 

e^,  Cp,  and  e^;  the  thermal  activation  energies,  - E^.  and  E.^,-  E^;  and  the 
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photoionization  rate  or  the  optical  emission  rates  of  electrons  and  holes,  e° 

and  e or  their  corresponding  photoionization  cross  sections,  0°  and  a°.  In 
P n P 

addition,  the  concentration  of  the  optically  sensitive  centers  directly  affects 

the  dtcector  sensitivity. 

The  detector  sensitivity  and  bandwidth  are  also  affected  by  the  pres- 
ence of  deep  level  recombination  centers  in  the  semiconductor  detector  material 
or  device  structures.  The  most  important  parameters  for  these  deep  levels  are 
the  thermal  capture  and  emission  rates,  thermal  activation  energies  and  concen- 
trations which  were  listed  above  for  the  optical  sensitive  centers. 

These  parameters  were  illustrated  in  the  energy  level  diagrams  in 
Figures  1.1,1  and  1.1.2. 

Although  many  experimental  studies  have  been  made  on  the  energy  level 
positions  of  the  imperfection  centers  in  silicon,  neither  precise  data  on  energy 
level  nor  on  the  kinetic  capture  and  emission  rate  parameters  indicated  above 
are  available.  These  parameters  are  essential  in  determining  the  optimum  detec- 
tor performance. 

The  purpose  of  the  studies  undertaken  in  this  project  was  to  apply  the 
transient  capacitance  techniques  to  obtain  the  precise  values  of  these  parameters 
for  potentially  useful  optical  centers  and  the  process  introduced  impurity  and 
defect  recombination  centers  under  a wide  range  of  temperature,  electric  fields, 
and  optical  spectra. 

The  transient  capacitance  techniques  have  been  reviewed  in  seven  invited 
lectures  given  by  C.  T.  Sah. [12-18]  The  texts  of  these  lectures  are  published 
in  four  papers  in  the  open  literature  with  emphasis  on  different  application 
aspects. [12,15,18,19] 

During  this  five-year  study,  important  and  powerful  extensions  of  the 
transient  capacitance  techniques  have  been  developed.  In  addition,  extensive 
applications  of  the  transient  capacitance  technique  were  also  made  to  study  the 
diffusion  and  generation-annealing  kinetics  of  the  deep  level  centers. 

3.1  EXPERIMENTAL  TECHNIQUES 

The  transient  capacitance  techniques,  first  extensively  used  by  Sah 
and  his  students  since  196")  [12],  have  been  extended  and  improved.  These  tech- 
niques include  [18,19]  the  VoZtagz  StinuZcutzcL  CapaCAXancz  (VSCAP),  previously 
known  as  the  dark  capacitance  transient;  the  TheAmotty  St^muZcuted  CapacuXancc 
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(TSCAP)  ; and  the  LLghjC  StXjnuZcutzd  CdpacAXajice.  {LSCA,P)  , previously  known  as  the 
photocapacitance.  These,  together  with  material  preparation  (growth  of  impurity- 
doped  crystal)  and  on-line  computer  control  of  experiments  and  analysis  of  experimen- 
tal data  are  described  in  this  section. 

3.1.1  EDGE  EFFECT  IN  TSCAP 

Electrons  and  holes  trapped  at  an  impurity  or  defect  center  in  the 
edge  region  of  the  junction  space  charge  region  can  give  rise  to  a change  in 
the  high-frequency  capacitance  of  the  junction  when  the  trapped  electrons  and  holes 
are  released  either  thermally  by  increasing  the  sample  temperature  or  optically 
with  a monochromatic  light.  The  effect  was  first  discovered  and  used  by  Walker 
and  Sah  in  1973  [20]  and  it  was  shown  that  one  can  observe  energy  levels  as 
shallow  as  60  meV  from  a band  edge  during  heating  from  77K  to  300K  in  a TSCAP 
experiment.  On  the  other  hand,  the  VSCAP  method  which  is  used  to  discharge  or 
charge  the  traps  in  the  space  charge  region,  can  only  reach  levels  no  shallower 
than  about  150  meV  from  one  of  the  two  band  edges  at  77K.  Thus,  the  edge  effect 
TSCAP  provides  a convenient  method  to  determine  the  presence  of  shallow  levels 
in  the  band  gap  without  having  to  use  a liquid  helium  dewar.  However,  temperatures 
lower  than  77K  is  still  needed  to  obtain  the  kinetic  parameters  of  very  shallow 
levels,  i.e.,  the  capture  and  emission  rates,  using  the  VSCAP  or  LSCAP  methods. 

The  theory  of  the  edge-effect  TSCAP  was  carefully  verified  experimen- 
tally by  Wang  [21];  in  particular,  the  dependences  of  the  magnitude  of  the 
capacitance  change  on  reverse  bias  voltage  was  verified.  Wang  also  applied  the 
method  to  determine  the  thermal  activation  energy  of  electrons  trapped  at  the 
Bi  level  in  silicon,  which  is  reported  in  the  subsection  on  Bi  in  this  chapter. 

3.1.2  EDGE  EFFECT  IN  VSCAP 

In  the  original  theory  of  the  voltage  stimulated  capacitance  (VSCAP) 
or  dark  transient  capacitance  of  a p/n  junction  described  in  our  1970  paper 
[22],  the  staircase  depletion  approximation  was  made  for  the  junction  space 
charge  region.  For  many  applications,  this  approximation  was  found  to  be  in- 
adequate because  of  the  contributions  from  the  edge  region.  Two  of  the  situa- 
tions in  which  the  edge  effect  becomes  important  are  investigated.  This  sub- 
section describes  one  of  these:  the  edge  effect  at  low  reverse  bias  in  profil- 
ing the  impurity  concentrations. 
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A detailed  theoretical  analysis  was  undertaken  by  Sah  and  Neugroschel 
In  1975  [23]  to  provide  a first  order  theory  of  the  contribution  of  the  edge 
region.  With  such  an  extension  of  the  theory.  Impurity  profiles  within  0.25 
micron  from  the  metal-semiconductor  Interface  of  a Schottky  barrier  or  the  junc- 
tion boundary  of  a p+/n  junction  can  be  obtained  with  good  reliability.  The 
Impurity  concentration  profiles  are  obtained  from  the  total  transient  capacitance 
change  with  reverse  bias  voltage.  In  order  to  get  the  concentration  profile 
near  the  junction  Interface,  VSCAP  data  at  low  reverse  or  even  small  forward  bias 
are  needed.  The  effect  of  the  edge  region  Is  large  at  low  reverse  bias  voltage 
and  it  significantly  reduces  the  computed  variation  of  the  impurity  concentration 
with  position  from  the  depletion  theory  which  completely  neglects  the  edge  region. 

The  first  order  theory  has  been  further  improved  by  Marchand  [24] , 
which  removed  the  approximation  of  constant  impurity  concentrations  In  the  edge 
region  used  by  Sah  and  Neugroschel  [23].  The  Improved  theory  is  now  used  rou- 
tinely in  a computer  program  to  evaluate  the  shallow  and  deep  level  Impurity 
concentration  profiles  In  silicon. 

Further  improvement  of  the  theory.  Including  the  Debye  length  in  the 
quasi-neutral  region,  will  be  made  to  complete  the  theoretical  development. 

3.1.3  CAPTURE  RATE  BY  VSCAP  IN  THE  EDGE  REGION 

The  effect  of  incomplete  filling  of  the  traps  which  were  subsequently 
emptied  to  give  the  VSCAP  was  first  observed  as  a filling  time  dependence  of 
the  total  VSCAP.  This  phenomenon  has  recently  been  analyzed  in  detail  and  the 
theory  verified  with  experimental  data  by  Marchand.  [24]  It  provides  a powerful 
new  method  of  measurement  of  the  thermal  capture  rate  of  electrons  and  holes  in 
the  edge  region  where  the  electric  field  is  low.  A proper  asymptotic  expansion 
of  the  theory  can  give  the  zero-field  or  near-thermal  equilibrium  capture  rates 
from  the  experimental  data.  The  thermal  equilibrium  capture  rates  are  the  para- 
meters necessary^  to  compute  the  carrier  recombination  lifetimes. 

3.1.4  SERIES  RESISTANCE  EFFECT  IN  VSCAP  MEASUREMENTS 

At  low  temperatures,  the  series  resistance  of  the  quasi-neutral  region 
and  the  edge  layer  of  a p/n  junction  or  Schottky  barrier  space  charge  region 
will  greatly  affect  the  transient  capacitance  measurements.  A detailed  inves- 
tigation of  the  series  resistance  effect  was  undertaken  by  Woodley  [25],  who 
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divided  the  edge  region  of  the  space  charge  layer  of  the  junction  into  several 
lumps.  Each  lump  is  represented  by  an  RC  parallel  equivalent  circuit  model. 

The  small-signal  capacitance  and  conductance  versus  temperature,  frequency  and 
junction  reverse  bias  voltage  are  obtained  and  compared  with  the  lump  circuit 
models.  Excellent  agreement  has  been  obtained  at  24K  as  indicated  in  Figure  3.1.1 
using  one  RC  lump  for  the  edge  region.  By  including  such  a model,  the  transient 
capacitance  data  at  very  low  temperatures  can  now  be  reduced  to  give  the  emis- 
sion and  capture  rate  constants  at  the  shallow  level  centers. 

3.1.5  SAMPLE  PREPARATION 

Two  unconventional  as  well  as  a number  of  conventional  device  fabri- 
cation techniques  were  developed  to  give  the  desired  samples  for  studying  the 
properties  of  shallow  and  deep  levels  in  silicon.  One  of  the  two  unconventional 
methods  was  the  development  of  crystal  growth  and  doping  techniques  to  provide 
silicon  crystals  which  are  doped  with  a shallow  and  a deep  level  impurity.  The 
doubly  doped  crystal  of  most  concern  to  this  study  was  the  silicon  crystal  doped 
with  both  phosphorus  and  bismuth  with  the  latter  at  a lower  concentration  to 
study  the  thermal  capture  and  emission  and  optical  emission  of  electrons  trapped 
at  the  B1  center.  The  second  development  was  the  successful  search  for  a metal 
which  would  provide  a good  rectifying  Schottky  barrier  on  p-type  silicon  from 
77  to  300K.  After  trials  with  many  impurities,  it  was  found  by  Tsao  that  Mg 
gives  a very  good  Schottky  barrier  even  at  room  temperature. [3]  Encapsulation 
or  storage  in  innert  ambient  was  found  to  be  necessary  to  prevent  deterioration 
of  the  Mg  on  p-Si  Schottky  barrier. 

3.1.6  ON-LINE  COMPUTER  CONTROL  AND  ANALYSIS  OF  EXPERIMENTS 

During  the  five-year  contract  period,  two  generations  of  minicomputer 
equipment  were  Installed  and  put  into  operation  for  on-line  control  and  analyses 
of  experiments.  The  on-line  computer  approach  is  necessary  to  relieve  the  point 
by  point  data  recording  and  analysis  methods,  which  are  tedious,  time-consuming, 
and  frequently  prone  to  human  error. 

The  first  generation  (acquisition  began  in  1970)  automatic  data  ac- 
quisition and  analysis  equipment  put  into  operation  in  our  laboratory  was  a 
programmable  desk  calculator  (Wang  720C)  with  a 500K  Byte  disc  and  a printer- 
plotter  terminal.  Part  of  the  A/D  and  D/A  interface  cards  as  well  as  the  entire 
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ireasuremeat  control  circuits  were  designed  and  built  by  T.  Y.  Lo  from  scratch, 
which  is  reported  in  his  doctoral  thesis.  After  successfully  installing  and 
operating  this  system,  it  became  evident  that  further  development  of  the  appli- 
cation programs  of  this  system  was  impractical  due  to  the  extremely  long  time 
(=  one  year  and  four  graduate  students)  required  to  write  and  debug  a program. 

This  drawback  came  from  the  low-level  calculator  language  used  by  the  Wang  720C 
computer  and  it  was  further  exemplified  by  the  limited  disk  storage  available. 

A second  generation  minicomputer— controlled  data  acquisition 
and  analyses  system  was. selected  to  provide  both  high  level  language  (FORTRAN  IV) , 
multiprogramming-multiuser  and  comprehensive  instrument-control  capabilities. 

A Hewlett-Packard  21MX  minicomputer  system  was  purchased  with  a 5M  Byte  disk, 
a 45  ips  tape  drive,  and  a Versatec  printer-plotter.  A number  of  instrumenta- 
tion Interfaces  have  been  or  are  being  added  to  control  digital  voltmeters  and 
d.c.  power  supplies.  An  Intel  8080  microcomputer  was  also  designed  and  built 
by  T.  Y.  Lo  to  serve  as  a stand-along  data  collection  terminal  but  programmed 
by  the  HP21MX  minicomputer  in  FORTRAN  IV.  This  microcomputer  will  enable  us  to 
free  the  minicomputer  for  other  tasks  while  a long-duration  experiment  is  being 
run  by  the  microcomputer. 


3.2 


EXPERIMENTAL  RESULTS 


Experimental  measurements  of  the  thermal  emission  and  capture  rates 
and  optical  emission  rates  are  described  for  phosphorus  centers,  bismuth  cen- 
ters, thermally  generated  centers,  impurity  centers  from  the  Boron+  source, 
gold  centers  and  a number  of  other  impurity  cneters  (Ag,  Co,  Cr,  S,  Ta  and  Zn) . 


3.2.1 
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The  phosphorus  energy  level  and  the  thermal  emission  rate  are  measured 
by  the  transient  capacitance  method  while  the  energy  level  Was  also  obtained 
from  conductance  measurements  of  a N+/N/N+  diode.  The  conductance  data  from  20 
to  33K  was  obtained  by  Nlng  and  Fu  (Figure  3.2.1)  and  fitted  to  the  Arrhenius 
equation  G = Aexp(-AE/kT)  and  to  give  a thermal  activation  energy  of  E*42  meV. 
The  difference  between  this  and  other  thermal  and  optical  energy  measurements 
which  gave  44  meV  was  probably  due  to  the  Inaccuracy  in  our  temperature  cali- 
bration and  measurements.  The  conductivity  measurement  is  a classical  method 
and  can  give  accurate  thermal  activation  energy  at  thermal  equilibrium  or  zero 
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Figure  3.2.2  Thermally  stimulated  capacitance  for  sulfur  doped  silicon 

diode.  Peaks  1 and  2 are  yet  unidentified.  Peak  3 is  due  to 
phosphorus  donors  in  the  edge  region  of  the  junction.  Peaks  4 
and  5 are  due  to  neutral  and  singly  ionized  sulfur  centers  in 
the  depletion  region  of  the  junction. 
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electric  field  when  experimental  errora  are  eliminated. 

The  phosphorus  energy  level  was  also  measured  by  the  TSCAP  method 
applied  to  a p+/n  junction  on  phosphorus  doped  silicon  substrate.  The  deriva- 
tive of  TSCAP  with  respect  to  temperature  is  shown  in  Figure  3.2.2.  The  major 
peak,  labeled  3,  at  42. 6K  comes  from  the  crossing  of  the  bulk  Fermi  level  and 
the  phosphorus  donor  ground  state  energy  at  44.5  meV  in  the  edge  region  of  the 
space  charge  layer  of  the  p+/n  junction.  The  smaller  and  broader  peaks  at 
higher  temperatures,  labeled  4 and  5,  are  from  the  two  sulfur  donor  levels  in 
the  upper  half  of  the  energy  gap  which  will  be  discussed  in  detail  under  Sulfur 
in  section  3.3.7.  The  other  two  smaller  fJeaks,  labeled  1 and  2,  at  34. 7K  and 
30. OK  have  not  been  identified  but  either  come  from  deionization  of  shallower 
levels  associated  with  unknown  impurities  or  from  the  excited  states  of  phos- 
phorus at  32  and  34  raeV  below  the  conduction  band  edge.  The  effect  of  the 
excited  state  on  the  dC/dT  peaks  was  analyzed  by  Woodley  in  which  the  edge  re- 
gion and  the  series  resistance  of  the  quasi-neutral  region  were  taken  into 
account.  Additional  dC/dT  measurements  were  also  made  on  diodes  with  different 
thermal  history,  which  showed  the  main  peak  but  a variety  of  smaller  peaks. 

Both  the  theoretical  and  experimental  results  suggest  that  the  smaller  peaks  at 
low  temperature  are  associated  with  the  other  Impurities  or  defects. 

The  thermal  emission  rate  and  activation  energy  of  electrons  trapped 
at  the  phosphorus  level  were  obtained  from  transient  capacitance  measurements 
using  a Al/n-Si  Schottky  barrier  diode  where  the  n-Si  was  doped  with  both  phos- 
phorus and  bismuth.  The  results  obtained  by  C.  T.  Wang  and  T.  J.  Woodley  are 
shown  in  Figure  3.2.3  and  fitted  to  the  Arrhenius  equation  e^  = Aexp(-AE/kT) , 
giving  A = (4±5)  x 10^  sec  ^ and  AE  = 53±5  meV.  This  experimental  thermal  ac- 
tivation energy,  approximately  53  meV,  is  higher  than  the  phosphorus  ground 
state  energy  of  44.5  meV.  The  discrepancy  is  probably  due  to  the  deionization 
effect  of  the  quasi-neutral  region  of  the  diode  resulting  in  a high  series  re- 
sistance which  distorts  the  capacitance  transient  considerably. 


3.2.2 


BISMUTH 


The  thermal  activation  energy  of  Bi  in  Si  was  first  measured  by  the 
VSCAP  method  using  a Al/n-Si  Schottky  barrier  diode  which  was  doubly  doped  with 
phosphorus  and  bismuth.  The  results  are  shown  in  Figure  3.2.3  giving  a pre- 
exponential  factor  of  l2t2)  x 10  sec  and  a thermal  activation  energy  of  67.1 
meV  compared  with  the  published  value  of  70.9  meV.  The  thermal  emission  rate 
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Figure  3.2.3  Reciprocal  time  constant  or  thermal  emission  rate  of  electrons 
trapped  at  the  phosphorus  and  bismuth  impurity  atoms  in  the 
depletion  region  of  the  reverse  biased  Schottky  Barrier  diode. 
Device  No.  101D2  pin  1,  V_=20  Volts. 


I 


of  crapped  electrons  obtained  in  this  experiment  was  not  available  previously. 

The  properties  of  Bi  were  also  obtained  by  the  VSCAP  both  in  the  de- 
pletion layer  and  the  edge  layer  of  the  space  charge  region  of  Al/n-Si  Schottky 
barrier  and  diffused  p+/n  diodes  on  silicon  doped  with  Bi  only.  The  position 
dependences  of  the  Bi  concentration  for  the  two  diodes  are  given  in  Figure  3.2.4. 
The  energy  level  of  bismuth  calculated  from  the  d(TSCAP)/dT  peak  due  to  ioniza- 
tion of  the  Bi  level  in  the  edge  region  was  74.5  meV  [21]  in  good  agreement  with 
optical  data  of  70.9  meV. 

The  Bi  level  is  also  obtained  from  the  temperature  dependence  of  the 
electron  concentration  computed  by  Woodley  from  his  capacitance  versus  tempera- 
ture data.  The  data  are  shown  in  Figure  3.2.5  which  compares  well  with  the 
theoretical  curve  of  a total  bismuth  concentration  of  = 3 x 10^**  Bi/cm^  and 
a Bi  ground  state  energy  level  of  = 70.90  meV.  The  experimental  data 

again  give  a somewhat  higher  activation  energy  than  70.90  meV. 

3.2.3  THERMALLY  GENERATED  CENTERS 

Deep  levels  introduced  inadvertently  during  high  temperature  process- 
ing of  silicon  devices  were  studied  in  detail  to  determine  the  origin  and  the 
source  of  these  centers.  A detailed  account  has  been  published  in  the  open  litera- 
ture [26]  as  well  as  in  Scientific  Report  No.  2.  In  summary,  there  are  four 
process  induced  or  thermally  generated  centers  which  are  donor-like  and  in  the 

upper  half  of  the  silicon  band  gap  and  which  have  high  concentrations  with  a 

13  -3 

U-shaped  concentration  profile  of  10  cm  in  the  flat  region  and  sharply  in- 
creasing values  within  a few  microns  of  the  two  surface  layers.  These  centers 
are  produced  when  the  sample  is  heated  above  1100®C  and  rapidly  cooled  to  room 
temperature  provided  that  there  is  a highly  disordered  surface  layer  on  the 
silicon  slice.  The  disordered  surface  layer  serves  as  the  source  of  these  re- 
combination centers.  It  can  be  either  a lapped  surface  or  a surface  with  a 
high  concentration  diffused  boron  or  phosphorus  layer.  Experimental  evidence 
showed  that  these  centers  are  vacancy  or  vacancy-impurity  complexes  and  that 
they  can  be  annealed  out  only  after  prolonged  and  slow  cooling  from  1100®C  to 
600“C.  A few  other  centers  were  also  observed  during  the  thermal  generation 
and  quenching  experiments  [26]  but  they  are  at  much  lower  concentrations  and 
are  generally  gettered  readily  by  the  n+  phosphorus  diffused  layer.  These 
levels  are  thought  to  be  associated  with  residual  impurities  introduced  during 
the  high  temperature  processes. 
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4 The  concentration  of  Bi  as  a function  of  distance  from  the  pn  junction 
which  is  0.5  micron  from  the  surface  for  the  diffused  p+/n  labeled  (2) 
and  for  the  Schottky  Barrier  diode  labeled  (1). 


.5  A comparison  of  the  electron  concentration  calculated  from  experi- 
mental data  of  high  frequency  SB  capacitance  versus  temperature  with 
the  theoretical  electron  concentration  of  silicon  doped  with 
3.0xl0i^  Bi/cm3  with  E„-E..=  70.90  meV. 


3.2.4 


IMPURITY  CENTERS  FROM  THE  BORON+  SOURCE 


When  the  Boron+  solid  source  is  used  to  produce  diffused  p+/n  silicon 

junctions  in  the  range  of  1000  to  1200C,  four  energy  levels  in  the  upper  half 

of  the  energy  gap  are  observed.  Detailed  VSCAP  measurements  of  the  thermal 

emission  rate  of  electrons  trapped  at  these  levels  are  obtained  by  R.  L.  Marchand 

and  A.  R.  Stivers  in  a sample  diffused  at  1200C.  The  data  are  fitted  to  the 

t 2 

Arrhenius  equation  e^  = A(T/300)  exp(-E/kT)  to  give  the  following  results; 


Table  2 

Parameters  of  Impurity  Centers  from  the  Boron+  Source 
-1  -3 

E(meV)  A(sec  ) N^(cm  ) Origin  (Tentative) 


93.  ± 2.2 

(7.7  ±2.1  )xl0‘* 

>6x10 

Mg° 

256.  ± 4. 

(7.6  ±2.8  )xl0ll 

6x10 12 

Mg"*" 

279.  ±13. 

(5.0  ±6.  )xl0ll 

3x10 12 

? 

534. 7±  1.3 

(1.04±0.07)x10^2 

2x10 13 

? 

A preliminary  search  of  the  lower  half  of  the  band  gap  gave  null  re- 
sults. The  two  shallowest  levels  given  in  the  above  table  appear  to  be  assoc- 
iated with  the  interstitial  Magnesium  donor  center  whose  energy  levels  were 

previously  determined.  The  detected  concentration  of  Mg*^  level,  N_  > 6x10^^, 

+ 

is  smaller  than  that  of  Mg  due  to  incomplete  Initial  filling  discussed  in 
section  3.1.3.  The  two  other  levels  are  probably  associated  with  the  impuri- 
ties present  in  the  Boron+  source.  From  the  data  sheet  of  the  Boron+  source, 
it  is  known  that  SrO,  MgO,  AI2O3  and  BaO  are  presented  as  impurities  at  a 

concentration  of  10%  or  more. 

At  a lower  diffusion  temperature  using  the  Boron+  source,  such  as 
llOOC,  the  impurity  concentrations,  remained  about  the  same.  However,  the 

279  meV  level  given  in  Table  2 above  is  missing.  Instead  a new  level  at  335t3 
meV  (^2x10^ ^cm"^)  was  detected. 

Detailed  studies  are  continuing  on  the  origin  and  generation-annealing- 
diffusion  kinetics  of  these  centers  as  well  as  the  feasibility  of  using  such  a 
multi- Impurity  diffusion  source  to  introduce  deep  and  shallow  level  impurity 
centers  with  potentially  useful  properties  in  detector  applications. 
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Figure  3.2,6  The  concentration  profile  of  the  aluminum  and  E^+392  meV 

hole  trap,  ^292'  phosphorus  concentration  profiles  are 

computed  from  expC-x^/AP-t)  using  three  slightly  different 
values  for  the  phosphorus  diffusivity  V. 


3.2.5 


ALUMINUM-OXYGEN  COMPLEXES 


Using  VSCAP,  a detailed  study  was  made  on  the  centers  thermally  gen- 
erated in  aluminun-doped  silicon  crystals,  grown  both  by  the  Czolchralski  and 
the  float  zone  methods.  Some  detailed  results  have  been  published  in  two  papers 
in  the  open  literature. [27,28]  In  summary,  three  energy  levels  in  the  lower 
half  of  the  silicon  energy  gap  and  one  energy  level  in  the  upper  half  of  the 
silicon  energy  gap  were  observed.  The  upper  level  was  found  to  be  located  near 
the  silicon  surface  while  the  three  lower  levels  were  distributed  deep  into  the 
silicon  bulk.  A detailed  study  of  the  diffusion  profile  indicates  that  these 
centers  are  associated  with  the  presence  of  oxygen  as  well  as  aluminum  in  the 
form  of  an  Al-0  pair  or  higher  aluminum-oxygen  complexes.  Experiments  at  1150C 
gave  an  out-diffusion  profile  of  the  deep  level  in  a sample  with  a high  phos- 
phorus-concentration diffused  surface  layer.  This  is  shown  in  Figure  3.2.6. 

The  out-dlf fusion  profile  of  the  392  meV  level  shown  in  the  inset  gave  a diffu- 
sivity  equal  to  that  observed  for  oxygen  in  silicon  at  1150C,  4.4x10”^ 'cm^/ sec. 
This  suggests  that  the  center  is  an  aluminum-oxygen  complex  which  dissociates 
at  1150C  and  the  oxygen  diffuses  into  the  high-concentration  phosphorus  diffused 
layer  which  serves  as  a sink  for  oxygen. 

The  oxygen-aluminum  complex  model  has  also  been  verified  by  experiments 
on  aluminum-doped  float-zone  silicon  crystals  which  have  low  oxygen  concentration. 

3.2.6  GOLD 

Gold  was  Che  most  commonly  used  impurity  in  silicon  to  reduce  the  min- 
ority carrier  lifetime  and  also  the  most  troublesome  impurity  that  can  be  easily 
introduced  into  the  silicon  device  during  high  temperature  processing. 

Detailed  thermal  emission  rate  data  of  the  electrons  and  holes  at  the 
two  energy  levels  have  been  published  in  the  open  literature  for  some  time  and  are 
reproduced  in  Figure  3. 2 . 7 . [18,19]  However,  the  thermal  capture  rates  had  not 
been  obtained  until  a preliminary  set  of  data  were  published  recently. [18,19] 

This  is  given  in  Figure  3.2.8.  Preliminary  capture  rate  data  has  also  been 
obtained  by  the  edge  effect  method  discussed  in  section  3.1.3.  This  is  given 
in  Figure  3.2.9.  One  of  the  main  problems  in  the  capture  rate  measurement  is 
that  the  observed  capture  process  occurs  in  a high  electric  field  region  where 
the  hot  electron  (or  hole)  effect  is  appreciable  while  the  capture  process  that 
controls  the  usual  electron-hole  recombination  lifetime  occurs  in  a low  field 
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second 


Figure  3.2 


.7  The  thermal  emission  rates  of  electrons,  e^, , and  holes,  e^„ 

nl  pO, 

trapped  at  the  gold  acceptor  center,  E^-547  mV  or  Ey+590  mV, 
and  the  thermal  emission  rates  of  holes,  ®p  trapped  at  the 
gold  donor  center,  Ey+345  mV,  obtained  from  the  dark  transient 
capacitance  measurements  as  well  as  the  dark  current  transient 
measurements,  the  latter  for  time  constants  less  than  1 second 


Figure  3.2.8  The  electric  field  and  temperature  dependences  of  the  thermal 

capture  coefficients  of  electrons  and  holes  at  the  gold  acceptor 
and  donor  centers  in  silicon  from  the  photo  capacitance 
experiments. 


or  quasi-neutral  region  which  is  essentially  at  thermal  equilibrium.  The  edge 
effect  method  for  measurement  of  the  capture  rate  discussed  in  section  3.1.3 
gives  essentially  the  thermal  equilibrium  capture  rate. 

The  concentration  profiles  of  gold  have  been  measured  within  a few 
micron  of  the  silicon  surface  of  the  Al-Si  interface  of  a Schottky  barrier.  A 
substantial  amount  of  gold  out-diffusion  into  the  aluminum  sink  has  been  dis- 
covered at  temperatures  as  low  as  450  C during  a 180  sec  die-bond  heating  of  the 
silicon  chip  on  a TO-18  header.  Illustrative  examples  are  shown  in  Figures  3.2.10 
and  3.2.11  [18,19].  In  the  silicon  sample  with  an  essentially  constant  initial 
gold  concentration  shown  in  Figure  3.2.10  U8],  the  gold  out-diffusion  into  a 
surface  sink  is  apparent.  In  Figure  3.2.11  [19],  two  phenomena  are  evident. 

First,  the  initial  gold  concentration  is  an  In-dif fusion  profile  which  is  ob- 
tained by  gold  diffusion  at  875C.  Second,  the  slight  decrease  of  the  gold  con- 
centration profile  near  the  surface  again  indicates  gold  out-diffusion  at  450C 
during  the  die-bonding  operation. 

3.2.7  DEFECT  CENTERS  FROM  ION  IMPLANTATION 

A large  number  of  electrically  active  lattice  defects  and  vacancy- 
impurity  complexes  are  formed  during  ion  implantation.  While  most  of  these 
centers  can  be  annealed  out  by  a subsequent  high  temperature  heating  above  about 
900C,  some  may  persist.  Extensive  measurements  of  the  thermal  emission  rates 
and  activation  energies  are  made  on  both  n-type  and  p-type  silicon  after  they 
are  exposed  to  energetic  ions,  including  lOOKeV  h"^,  ISOKeV  300KeV  2®Si^, 

1.3MeV  ®'*S  , and  2MeV  ^®0^  for  a total  flux  of  10^^  to  10^®  atoms/cm  . Prior  to 
the  ion  implantation,  very  shallow  p+/n  and  n+/p  diffused  diodes  were  made  with 
appropriate  diffused  and  MOS  guard  rings  to  cut  off  surface  channels  and  peri- 
meter breakdown  effects.  Some  diodes  are  heat  treated  in  a dry  nitrogen  ambient 
prior  to  the  transient  capacitance  measurements. 

Typical  results  are  shown  in  Figure  3.2.12  for  samples  Implanted  with 
10^^  0/cm^.  The  labels  are  as  follows.  The  labels  0,  01,  and  04,  at  the  top 
of  the  curve.  Indicate  a subsequent  anneal  of  90  seconds  at  150*C  for  01,  400*C 
for  04  and  no  annealing  for  label  0.  The  numerical  labels  on  the  side  of  these 
curves  are  the  peak  temperatures  of  the  d(TSCAP)/dT  curve  and  the  thermal  acti- 
vation energies  relative  to  the  indicated  band  edge  using  m>0  in  the  Arrhenius 
equation  given  by  e * Aexp(-E/kT).  The  concentration  of  these  centers  is  also 
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indicated  on  the  side  o£  the  lines  but  they  are  some  average  values  since  ion 
generated  GRT  centers  are  distributed  very  nonunlfomly  into  the  silicon  sur- 
face. Detailed  measurements  of  the  GRT  concentration  profiles  are  incomplete 
at  this  time.  For  comparison,  several  centers  generated  by  one  MeV  electron 
irradiation  are  also  given  in  Figure  3.2.12  with  a label  'e',  showing  similari- 
ties between  them  and  the  ion  generated  centers.  In  addition,  a center  gener- 
ated by  Si  ion  with  =*  170K  is  also  shown  which  coincides  with  one  of  the 

n 

centers  generated  by  the  one  MeV  electron. 

The  centers  observed  after  0 ion  Implant  are  generally  also  observed 
after  implantation  by  the  other  ions  just  mentioned.  This  result  indicates 
that  most  of  the  ion  generated  GRT  centers  are  not  related  to  the  ion  species 
used  but  are  lattice  defects,  vacancy  clusters  and  vacancy-impurity  complexes. 

It  is  evident  from  the  preliminary  results  given  in  Figure  3.2.12  that  there  is 
much  interesting  and  very  application-relevant  work,  such  as  integrated  IR 
detector  arrays,  that  needs  to  be  done  by  detailed  measurements  of  the  thermal 
activation  energy  and  emission  rates  as  well  as  the  concentration  profiles  of 
the  centers  generated  by  energetic  ions  during  ion  implantation  processes. 

3.2.8  OTHER  IMPURITY  CENTERS 

A number  of  other  impurity  centers  were  investigated,  including  Ag, 

Co,  Cr,  S,  Ta,  and  Zn.  The  results  are  summarized  below  since  all  are  in  open 
literature  except  Cr,  whose  preliminary  results  will  be  reported  here.  In  this 
summary,  the  thermal  emission  rates  are  fitted  to  the  Arrhenius  equation, 
e = A(T/300)^exp(-E/kT)  where  the  thermal  activation  energy  and  the  pre-exponen- 
tial factor.  A,  are  quoted  in  the  following  discussion.  Photolonlzatlon  proper- 
ties are  also  summarized  when  available. 

Silver 

The  thermal  emission  rate  measurements  of  the  Ag  impurity  in  silicon 
showed  that  it  is  amphoteric  with  a donor  level  located  at  378  meV  above  the 
valence  band  edge  and  an  acceptor  level  located  at  544  meV  below  the  conduction 
band  edge. [29]  Doping  of  silicon  by  Ag  diffusion  is  messy  because  Ag  tends  to 
flow  over  the  entire  Si  surfaces  and  alloys  with  Si.  The  thermal  emission  rates 
of  the  two  centers  are  given  below: 

Donor  Center:  E.^  - E^  » 378  ± 2 meV  A * (2.60±0.27)  xlO^^  sec~^ 

Acceptor  Center:  - E.j,  « 544  ± 2 meV  A = (0.31±0.026)xl0^  ^ sec~' 
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Figure  3.2.12  The  reciprocal  emission  rate  or  Che  transient  capacitance  decay 

time  constants  of  the  energy  levels  generated  during  Ion  implanta- 
tion  as  a function  of  reciprocal  temperature.-  The  data  are  fitted 
to  Che  Arrhenius  equation  e*Aexp(-E/kT) . 


ionized  sulfur  centers. 
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j Since  the  two  silver  energy  levels  are  very  similar  to  those 

i of  gold,  careful  measurements  were  made  to  compare  the  emission  rate  results 

I [30]  and  it  was  demonstrated  that  indeed  the  thermal  emission  rates  at  the  two 

silver  centers  can  be  distinguished  from  those  of  the  two  gold  centers. 

Cobalt 

Thermal  emission  rate  measurements  made  on  Co  diffused  silicon  showed 
that  Co  is  amphoteric  with  a donor  level  located  at  352  meV  above  the  valence 
band  edge  and  an  acceptor  level  located  at  534  meV  below  the  conduction  band 
edge.  The  thermal  emission  rates  of  the  two  centers  are  given  below: [30,31] 

Donor  Center:  = 352  ± 1 meV  A » (2.46±0.19)xl0^^  sec”^ 

Acceptor  Center:  E^,  - E^  + 534  ± 3 meV  A » (0.77±0.10)xl0^^  sec”^ 

A careful  check  against  the  thermal  emission  rates  of  the  gold  centers 
was  also  made  to  ascertain  that  the  two  Co  centers  were  not  due  to  gold  contam- 
inations. [30] 

Chromium 

Chromium  was  Introduced  into  silicon  during  Czochralskl  crystal  growth 
with  a segregation  coefficient  of  about  10"^  at  the  silicon  melting  point.  From 
thermal  emission  rate  measurements,  it  was  found  that  Cr  appears  to  Introduce  two 
energy  levels  into  silicon.  Preliminary  data  indicates  that  there  is  a shallow 
level  about  180  meV  from  the  valence  band  edge  and  a deep  level  at  232  meV  from 
the  conduction  band  edge.  Both  levels  appear  to  be  donor-like.  Accurate  meas- 
urements were  made  for  the  upper  donor  level  and  the  thermal  emission  rate  data 
is  given  by: 

Donor  Center:  E^  - E^  = 232  ± 3 meV  A * 7.6  x 10^^  sec~^ 

Additional  measurements  are  planned  to  obtain  more  accurate  data  on 
the  Cr  centers  in  silicon. 

Sulfur 

Very  detailed  studies  had  been  made  prior  to  the  start  of  this  five- 
year  contract  and  continued  during  this  contract.  Both  thermal  emission  rate 
and  photolonlzatlon  rate  measurements  were  made  as  a function  of  temperature, 
electric  field  and  photon  energy.  Two  donor  levels  were  found  which  are  consis- 
tent with  the  valence  bond  model  since  sulfur  has  six  valence  electrons  and  only 
four  are  needed  to  satisfy  the  silicon  covalent  bonds. 

The  thermal  activation  energy  of  the  donor  electrons  at  the  two  S cen- 
ters are  given  below. [32] 
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First  Donor:  = 276  ± 1 meV  A - (l,64±0.12)xl0*®  sec”^ 

Second  Donor:  E^  - E^  = 528  ± 4 meV  A = (1.03±0. 22)xl0^^  sec"^ 

Electric  field  dependences  of  the  thermal  capture  rates  of  electrons 
at  both  of  the  sulfur  donor  centers  were  measured  and  these  are  given  in  Figures 
3.2.13  and  3.2.14.  In  addition,  impact  ionization  rates  of  the  neutral  sulfur 
center  by  electrons  was  also  measured  and  shown  in  Figure  3.2.15. 

The  fine  details  of  the  photoionization  rates  were  first  obtained  for 
the  second  donor  by  Rosier  [33,34],  Indicating  a strong  photothermal  ionization 
process  below  the  ground  state  threshold  of  613  meV,  a value  obtained  from  ab- 
sorption niea3ure;i]ents.  The  photothermal  ionization  process  is  a two  step  pro- 
cess in  which  the  electron  is  first  optically  excited  from  the  613  meV  ground 
state  to  an  excited  state  and  then  thermally  released  from  the  excited  state 
into  the  conduction  band.  Fine  structures  were  observed  as  follows:  [see 
Flg-ire  3.3.12  In  section  3.4.3  for  the  experimental  data] 

524  meV  Is(Ai)  ls(T2) 

570  meV  Is(A^)  -*•  2po 

574  meV  IsCA^)  2s  (Tj) 

591  meV  Is(A^)  -►  2p^ 

where  the  ground  state,  Is(A^),  is  located  at  613  meV  below  the  conduction  band 
edge. 

Detailed  photolonlzatlon  rate  measurements  were  also  carried  out  for 
the  first  donor  center  which  contains  two  electrons. [35]  The  low  temperature  threshold 
obtained  by  extrapolating  the  data  of  59. 8K.  appears  to  reach  280  meV,  which  is 
nearly  chat  of  Che  thermal  activation  energy  from  thermal  emission  rate  meas- 
urements at  high  temperatures,  276  meV.  An  absorption  peak  was  observed  at 
285,4  meV  which  is  attributed  to  the  [Is(Aj),  Is(Aj)]  -*■  [IsCA^),  Isd^)]  optical 
transition  while  the  continuum  is  315  meV  above  the  ground  state  [Is (Aj ) , Is(Aj)]. 

The  broadening  of  the  absorption  line  at  285.4  meV  was  shown  to  be  due  to  lat- 
tice or  phonon  scattering  and  follow  a Lorentzlan  line  shape  (See  section  3.3.3). 

Tantalum 

Thermal  emission  measurements  on  Ta-doped  silicon  crystals  indicated 
Chat  there  are  two  donor  levels  below  Che  silicon  conduction  band  edge  with 
the  following  parameters: 
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; <E>  (lo'*  volts/cm) 

Figure  3.2.14  Field  dependence  of  Che  capture  rate  of  electrons  by  doubly 
Ionized  sulfur  centers. 


( 


Figure  3<.2.15  Field  dependence  of  the  impact  ionization  race  of  neutral 
sulfur  centers. 


-44- 


-A  3 


First  Donor:  E^,  - = 232  ± 5 meV  A = (1.50±0.83)xl0^®  sec“* 

Second  Donor:  E^  - E^  = 472  ± 8 raeV  A (1.47±0.78)xl0^®  sec"^ 

There  were  serious  difficulties  in  getting  sufficient  concentration  of  Ta  into 

Si  during  crystal  growth  due  to  the  very  low  segregation  coefficient,  estimated 
at  10“ ® or  less. 

Zinc 

Detailed  measurements  of  the  thermal  emission  races,  thermal  capture 
rates  and  photolonlzaClon  cross  sections  were  made  on  Zn-dlf fused  silicon  n+/p 
and  p+/n  junctions  [36,37,38]  because  of  the  considerable  interest  shown  by 
previous  workers  in  Zinc  as  a possible  IR  detection  center  in  silicon. 

Thermal  emission  rate  measurements  indicated  that  zinc  gives  two  accep- 
tor levels  in  the  silicon  band  gap,  located  at  297  and  571  meV  above  the  valence 
band  edge.  These  are  consistent  with  the  valence  bond  model  since  zinc  has  only 
two  valence  electrons  and  needs  two  additional  electrons  to  complete  the  tetrahedral 
covalent  bond  of  the  silicon  lattice.  A shallow  level  at  130  meV  from  the  val-  i 

ence  band  edge  is  also  observed  in  zinc  diffused  and  boron-doped  silicon  crys-  | 

) 

tal  which  was  shown  to  be  a BZna  complex.  The  thermal  emission  parameters  are  j 

listed  below. 

First  Acceptor:  E^  - E^  =■  297  ± 1 meV  A = (6.78±0.68  )xl0*^  sec“^ 

Second  Acceptor:  E^  - E^  = 571  ± 1 meV  A = (0. 14±0.002)xl0^^  sec~^ 

BZn3  Acceptor:  E^  - E^  ■ 167  ilOmeV  A = (1.5  ±2.  )xl0^°  sec~^ 

The  thermal  capture  rates  measured  in  high  electric  field  and  the 
electron  Impact  ionization  rate  are  also  measured  and  these  are  given  in 
Figures  3.2.16,  3.2.17,  and  3.2.18. 

. 

The  photoionization  cross  sections  show  some  temperature  and  field  de- 
pendences below  the  continuum  threshold  but  no  distinct  structures  or  peaks 
such  as  those  observed  for  the  two  sulfur  donor  levels. 

3.3  THEORETICAL  ANALYSES  AND  INTERPRETATIONS 

Theoretical  analyses  of  the  photoionization  cross  sections  of  Impurity 
centers  were  carried  out  for  the  following  cases:  the  Coulomb  and  the  square 
three-dimensional  impurity  potential  wells  with  one  bound  electron  (or  hole) 
and  the  Coulomb  potential  well  with  two  bound  electrons.  The  results  are  sum- 
marized In  the  following  subsections  and  compared  with  some  available  experi- 
mental data. 
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Figure  3.2.16  Electric  field  dependence  of  the  thermal  Figure  3.2.17  Electric  field  dependence  of  the  thermal 

capture  rate  of  holes  at  singly  Ionized  capture  rate  of  holes  at  doubly  ionized 

zinc  centers.  zinc  centers  at  the  temperatures  Indicated 
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Figure  3.2.18  Temperature  dependence  of  the  thermal  capture  rate  of 

electrons  at  singly  ionized  zinc  centers  and  Auger-impact 
emission  rate  of  electrons  at  doubly  ionized  zinc  centers 
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3.3.1  COULOMB  WELL  FOR  SHALLOW  LEVEL  CENTERS 

The  bound-state  wavefunctlon  is  given  by 

|i>  » (wa^)  exp(-r/a)  (3.1) 

where  a is  the  Bohr  radius.  The  conduction  band  state  is  given  by  the  plane 
wave 

lf>  *15  exp(ik*r)  (3.2) 

where  U is  the  volume  of  normalization.  The  photo ionization  cross-section  is 
calculated  from 

o(hm)  - t(E^ff/E^)2(n/e^)](4Tr2/3)a)4a)J^l<f|?li>|26(E^+Ej.-lla))  (3.3) 

where 

E -^/E  » effective  electric  field  ratio  at  the  impurity  ion  center 

eff  o 

n * optical  index  of  refraction 
* optical  dielectric  constant 
a * e^/Hc,  the  fine  structure  constant 
Ej  = ionization  energy  of  the  bound  state 
Using  the  initial  and  final  state  wave  functions.  Equations  (3.1)  and  (3.2), 
to  calculate  the  dipole  matrix  element  in  Equation  (3.3),  the  photoionization 
cross-section  is  given  by 

o(ha))  - [(E  ,,/E  )2(n/e  ) ] (4Tra/3)a2s(hw)  (3.4) 

et  t o o 

Here  the  dimensionless  shape  function  is  given  by 

S(X)  • 256AX[A(X-Ej,)]^^^[l+A(X-Ej.)]"^  (3.5) 

where  X = Hw,  A » 2m*a^/K^  and  the  spherical  band  E = H^k^/2m*  ■ Km  - E^  is 
used.  This  was  obtained  previously. [39] 

Equations  (3.4)  and  (3.5)  are  fitted  to  the  published  photoionization 
spectra  of  group  III  acceptors  in  silicon  using  the'  Bohr  radius  'a'  as  the  ad- 
justable parameter.  These  are  shown  in  Figures  3.3.1  to  3.3.4. 

The  effective  radii  'a'  thus  determined  are  plotted  in  Figure  3.3.5 
as  a function  of  E^.  In  Figure  3.3.5,  the  'a'  from  variational  calculations 
of  the  group  V donors  in  silicon  are  also  plotted.  All  the  'a'  values  thus 
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gure  3.3.3  Spectral  dependence  of  the  photolo,nlza-  Figure  3.3.4  Spectral  dependence  of  the  photoloniza' 

tion  cross-section  of  Ga  In  SI.  The  tlon  cross-section  of  In  In  SI.  The 

theory  Is  based  on  a Coulomb  potential  theory  Is  based  on  a Coulomb  potential 

well.  well. 


obtained  ace  on  the  line  E^a  > constant  as  required  by  the  Coulomb  model. 

3.3.2  SQUARE  WELL  FOR  DEEP  LEVEL  CENTERS 

Since  the  delta-function  potential  well  has  been  generally  used  to  fit 
the  observed  photoionization  cross-section  spectra  of  deep  level  impurity  with 
only  moderate  success,  this  delta  function  well  model  is  extended  to  a well  of 
finite  radius  and  well  depth.  The  square  well  is  a particularly  good  approxi- 
mation for  neutral  defect  centers  such  as  the  one-MeV  electron  irradiation  gen- 
erated centers  in  silicon. [40] 

The  spherical  square-well  potential  is  given  by 


V(r)  - -V, 


V(r)  - 0 


. (3.6) 


and  the  wave  function  for  a bound  S state  (i~0)  is 


|i>  - Ro(r)  Yoo(r) 


Ajo(r/e+vo)  Yoo(r)  r<a 

Bh^^^dr/T”)  Yoo(r) 


(3.7) 


where  a is  the  radius  of  the  well,  Jq  is  the  spherical  Bessel  function,  hg  ' is 
the  spherical  Hankel  function  of  the  first  kind,  and  Yqo  is  the  spherical  har- 
monic. e > 2m*E/h^  is  the  normalized  energy  and  vg  ~ 2m*Vo/h^  is  normalized 
well  depth.  The  constants  A and  B are  determined  by  normalization  and  conti- 
nuity of  the  wave  function  |l>  at  r>a.  The  final  states  are  taken  to  be  plane 
wave  states  given  by  Equation  (3.2)  which  can  be  expanded  in  spherical  harmonics 


|f>  » U ^^^exp(i^*r) 

- y’^^^47r^Il^jj^(kr)Y*^(r)^(^) 

£m 

to  facilitate  the  calculation  of  the  dipole  matrix  element  contained  in 
Equation  (3.3)  which  is 


(3.8) 


|<f|r|l>|2  - (Arr/U)  r3jj(kr)  Ro(r)dr 


(3.9) 
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The  final  result  of  the  photoionization  cross-section  after  integrating  over 
all  k of  the  final  states  in  Equation  (3.3)  is  given  again  by  Equation  (3.4) 
and  the  shape  function  is  now 


S(Hu))  -(2m*)lu»k/)l2a2) 


r3j^(kr)Ro(r)dr 

0 


(3.10) 


where  k^  * (2in*/ll^)  (Ku-Ej) . In  the  limit  of  ka-K). 

S(Jlu))  - (4H2/m*a2)Ej^^(Haj-Ej)^^^/(hw)3  (3.11) 

and  the  photoionization  cross-section  becomes 


a(Ku)  - aj^(Ku) 

- [(E^^j/E^)2(n/e^)](16Trah2/3m*)Ej^^(Ho»-Ej)^^^/(Huj)3  (3.12) 

which  is  the  familiar  result  of  the  delta-function  well  given  by  Lucovsky. [17] 
The  Integral  of  Equation  (3.9)  can  be  evaluated  in  close  form  but  is 
long  and  complicated.  A family  of  normalized  curves  is  shown  in  Figure  3.3.6 
to  Illustrate  the  shape  of  the  photoionization  spectra  as  a function  of  well 
size.  It  is  evident  that  the  spectra  becomes  more  sharply  peaked  and  the  peak 
moves  toward  the  threshold  as  the  well  radius  is  Increased.  Such  a sharply 
peaked  spectra  has  been  observed  for  some  one-MeV  electron  irradiation  Induced 
centers  in  silicon  as  indicated  in  Figure  3.3.7.  This  is  consistent  with  the 
model  that  these  electron  irradiation  Induced  centers  are  extended  lattice  de- 
fects which  are  electrically  neutral  and  can  be  characterized  by  a square  well. 
[40] 

In  Figure  3.3.8,  the  photolonlzatlon  data  of  In  in  silicon  is  fitted 
to  a square  well.  This  is  as  good  a fit  as  the  Coulomb  well,  shown  in 
Figure  3.3.4  and  both  the  square  well  and  the  Coulomb  well  fit  much  better 
than  the  delta- function  well.  This  result  Illustrates  that  the  photoionization 
spectra  Is  mainly  Influenced  by  the  bound  state  wave  function  away  from  the 
imperfection  center  and  it  does  not  probe  the  central  cell  part  of  the  imper- 
fection potential  and  cannot  distinguish  a Coulomb  well  from  a square  well. 

It  does  show  that  the  commonly  used  delta-function  well  is  oversimplified  and 
too  compact  for  most  imperfection  centers. 
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3.3.3 


C0UL(»1B  WELL  FOR  TWO-ELECTRON  CENTERS 


I 


The  photoionization  cross-section  for  a two-electron  center  was  ob- 
tained using  the  Helium-like  model. [35]  It  was  shown  that  the  experimental 
data  exhibit  a Lorentzlan  line  shape. 

In  the  electric  dipole  approximation  the  photoionization  cross-section 
for  a two  electron  neutral  center  is  given  by  an  equation  similar  to  (3.3): 

- [(Egjj/E^)2(n/e^)](4^2a/3)Hu)|l<f|?^+r^|i>126(Ej-E^-ttu)  (3.13) 

where  n is  the  optical  index  of  refraction,  is  the  optical  dielectric  con- 

stant, E ../E  is  the  effective  electric  field  at  the  center,  a=e^/lic  is  the 
ett  o 

fine  structure  constant.  Here  |i>  is  the  initial  state  or  ground  state  of  the 
neutral  center,  and  |f>  is  the  final  state  where  one  of  the  electrons  is  in 
the  continuum  state  and  the  other  electron  is  in  the  ground  state  of  the  singly 
ionized  center.  In  the  helium-like  model,  the  initial  ground  state  for  the 
two-electron  center  is  assumed  to  be 

|i>  - (l/ira^)  exp[r(rj+r2)/a]  (3.14) 


and  the  final  continuum  state  is  written  as 


|f>  » |k>  * (2Trb^)  ^^^{exp[(-ri/b)+ik*r2] 

+ exp [ (-i2/b)<-i^*r  j ] } (3.15) 

where  a is  the  effective  Bohr  radius  of  the  two-electron  neutral  center,  and 
b is  the  radius  of  the  one-electron  singly  ionized  center.  ^ Is  the  volume  of 
normalization  and  k is  the  wave  vector  which  characterizes  the  electron  in  the 
continuum.  The  summation  in  Equation  (3.13)  then  becomes  a sum  over  a].l  the  k 
values.  Using  Equations  (3.14)  and  (3.15)  to  calculate  the  electric  dipole 
matrix  elements  in  Equation  (3.13),  then 

aJ°^(h(o)  - (E^j^/E^)2(n/e^)(47ra/3)a2s^°^()iio)  (3.16a) 


where  the  spectral  dependence  is  contained  in  the  dimensionless  shape  function 

]3 


512 


2m*a2 

~Y? 


(b/a+1) 


6 [(2m*a2/Jl2)(j(u,-E5°^]^^^ 


[l+(2m*a2/lll2)(hui-E;^^)] 


(0),,6 


(3.16b) 


Lorenizion 


Here  is  the  one-electron  ionization  energy  of  the  neutral  center. 

After  the  ionization  of  one  of  the  electrons  from  the  helium-like 
neutral  center,  the  center  becomes  singly  ionized  and  is  hydrogen-like.  The 
photoionization  cross-section  for  the  singly  ionized  or  one-electron  center 
can  be  calculated  from  Equation  (3.3)  using  1 i>  = /l/irb^  exp(-r/b)  for  the 
initial  ground  state,  and  |f>  = lk>  = /1/U  exp(4ik*r)  for  the  final  continuum 
state.  The  results  were  given  by  Equations  (3.4)  and  (3.5)  and  with  proper 
notation  changes  we  have 

with 


=eff 

2 

11 

E 

e 

o 

Atto 


(yiu)) 


(3.17a) 


(Mm) 


256 


2m*b^ 

h^ 


[(2m*b2/H2)(Hm_E(‘^^)]3/2 

[ l+(2m*b2  /Ji2  ) (Idw-E^'*’^  ) ] ^ 


(3.17b) 


where  E^^^  is  the  ionization  energy  of  the  singly  ionized  center. 

In  order  to  evaluate  the  cross-sections  from  Equations  (3.16)  and 
(3.17)  for  a given  center,  we  need  to  know  the  effective  Bohr  radii  a and  b 
for  the  neutral  and  the  singly  ionized  centers.  At  present  these  values  are 
available  only  for  sulfur  centers  in  silicon. [42]  From  variational  calcula- 
tions of  the  energy  levels,  the  effective  Bohr  radii  for  sulfur  centers  in 
silicon  were  a ^ 8.1  a.u.  and  b ==  7.3  a.u. 

A comparison  of  theory  with  experiment  is  given  in  Figure  3.3.9  for 
the  photoionization  cross  section  of  the  two-electron  sulfur  center.  The  ab- 
sorption tail  near  and  below  the  threshold  energy,  315  meV,  comes  from  the  two- 
step  photothermal  ionization  processes  previously  discussed  in  section  3.3.7 
under  Sulfur,  which  has  a Lorentzlan  absorption  llneshape  as  indicated  in 
Figure  3.3.10. 

A comparison  of  the  theory  with  experiment  is  also  made  in  Figure  3.3.11 
for  the  photoionization  cross  section  of  the  one-electron  sulfur  center.  Again 
the  absorption  tall  near  and  belcw  the  threshold  energy,  613  meV,  comes  from 
the  many  two-step  photothermal  ionization  processes  previously  discussed  in 
section  3.3.7  under  Sulfur.  The  details  of  the  absorption  spectra  near  and  be- 
low the  threshold  are  shown  in  Figure  3.3.12. 
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The  material  parameters  which  control  the  performance  of  an  extrinsic 
semiconductor  IR  detector  are  delineated.  These  are:  (i)  the  carrier  mobili- 
ties, (ii)  the  thermal  capture  and  emission  rates  and  the  optical  emission 
rates  or  photolonlzatlon  cross-sections  of  electrons  and  holes  trapped  at  the 
recombination  and  optical  centers,  and  (iii)  the  concentration  of  these  cen- 
ters. These  centers  are  either  impurities,  defects,  or  impurity-defect  com- 
plexes. 

Mobility  studies  described  in  chapter  II  of  this  report  show  that 
electron  and  hole  mobilities  at  high  impurity  concentrations  are  limited  by 
scattering  of  the  randomly  located  impurity  ions.  A comparison  of  the  experi- 
mental data  with  several  theoretical  models  indicated  that  there  is  no  ion 
position  correlation.  Position  correlation  would  reduce  random  scattering  and 
Increase  mobility,  allowing  the  use  of  high  concentrations  of  optically  active 
Ions  to  increase  IR  detector  sensitivity. 

Fabrication  of  silicon  extrinsic  IR  detectors  requires  high  tempera- 
ture processing  during  which  impurities  and  defects  are  inevitably  Introduced. 
The  presence  of  minute  amounts  of  impurities  and  defects  will  drastically  affect 
the  detector  characteristics  such  as  detection  sensitivity,  bandwidth  or  response 
speed,  radiation  sensitivity  and  life.  Previous  classical  techniques,  such  as 
the  Hall  effect  and  the  photoconductivity  methods,  are  not  sensitive  nor  unam- 
biguous enough  to  detect  the  presence  of  these  impurity  and  defect  centers. 

New  methods,  using  the  transient  capacitance  measurement  of  a semiconductor  space 
charge  layer,  have  been  developed  to  detect  and  characterize  these  new  centers. 


A summary  of  these  methods  is  given  in  section  3.1. 

Many  impurity  and  defect  centers  in  silicon  have  been  studled'and  re- 
ported in  section  3.2,  including  phosphorus,  bismuth,  gold,  magnesium,  silver, 
cobalt,  chromium,  sulfur,  tantalum,  zinc,  aluminum-oxygen  complexes,  thermally 
generated  donors,  and  defect  centers  generated  by  energetic  ions.  In  principle, 
the  impurity  centers  can  be  introduced  or  eliminated  in  controlled  ways  during 
device  fabrication.  However,  the  impurity-oxygen  complexes  and  the  centers 
generated  thermally  during  high  temperature  fabrication  processes  and  the  cen- 
ters generated  during  ion  implantation,  if  not  eliminated  or  controlled,  may 
become  detrimental  to  the  silicon  extrinsic  IR  detector  performance.  Thus,  it 
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seems  highly  desirable  to  undertake  a systematic  measurement  program  to  charac- 
terize these  impurity  and  defect  centers  using  the  highly  sensitive  capacitance 
transient  techniques.  Such  a study  may  provide  not  only  information  on  new 
optical  centers  useful  in  extrinsic  silicon  detector  applications,  but  also  indi 
cate  fabrication  conditions  to  make  their  introduction  and  elimination  predict- 
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METRIC  SYSTEM 


BASE  UNITS: 

Quantity 


SI  PREFV/.S: 


Multiplication  Factors 


I 000  000  000  000 
1 000  000  000 
I 000  000 
1 000 
100 
10 
0.1 
0.01 
0.001 
0.000  001 
0.000  000  001 
0.000  000  000  001 
0.000  000  000  000  001 
0.000  000  000  000  000  001 


’ To  tw  avoidad  when  poaaible. 


SI  Symbol 


Formula 


langth 

metre 

m 

maaa 

kilogram 

kg 

tima 

second 

s 

alactric  cuirant 

ampere 

A 

thannodynamic  temparatuia 

kelvin 

K 

amount  of  substance 

mole 

mol 

luminous  intanaity 

candela 

cd 

SUPPLEMENTARY  UNITS: 

plana  angla 

radian 

rad 

.4* 

solid  angle 

steradian 

sr 

... 

DERIVED  UNITS: 

Accaiantion 

metre  per  second  squared 

m/s 

activity  (of  a radioactive  source) 

disintegration  per  second 

(disintegration|/s 

angular  acceleration 

radian  per  second  squared 

rad/s 

angular  velocity 

radian  per  second 

rad/s 

aree 

square  metre 

m 

density 

kilogram  per  cubic  metre 

kg/m 

electric  capacitance 

farad 

F 

A-sA/ 

electrical  conductance 

siemens 

S 

A/V 

electric  field  strength 

volt  per  metre 

V/m 

electric  inductance 

henry 

H 

V-s/A 

electric  potential  difference 

volt 

V 

W/A 

electric  resistance 

ohm 

V(A 

electromotive  force 

volt 

V 

W/A 

energy 

joule 

1 

N*m 

entropy 

joule  per  kelvin 

... 

l/K 

force 

newton 

N 

kg-m/s 

ftequency 

hertz 

Hz 

(cycle)/8 

illuminance 

lux 

lx 

im/m 

luminance 

candela  per  square  metre 

... 

cd/m 

luminous  flux 

lumen 

Im 

cd'sr 

magnetic  field  strength 

ampere  per  metre 

... 

A/m 

magnetic  flux 

we^ 

Wb 

V-s 

magnetic  flux  density 

tesla 

T 

Wb/m 

magnetomotive  force 

ampere 

A 

i/i 

power 

watt 

W 

pressure 

pascal 

Pa 

N/m 

quantity  of  electricity 

coulomb 

C 

A-s 

quantity  of  heat 

joule 

1 

N*m 

radiant  intensity 

watt  per  steradian 

W/sr 

specific  heat 

joule  per  kilogram-kelvin 

)/kg.K 

stress 

pascal 

Pa 

N/m 

thermal  conductivity 

watt  per  metre-kelvin 

W/m-K 

velocity 

metre  per  second 

m/s 

viscosity,  dynamic 

pascal-second 

... 

Pa-s 

viscosity,  kinematic 

square  metre  per  second 

m/s 

voltage 

volt 

V 

W/A 

volume 

cubic  metre 

m 

wavenumber 

reciprocal  metre 

(wave)/m 

work 

joule 

r 

N-m 

.SI  Symbol 
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